NORMALLY ELLIPTIC SINGULAR PERTURBATIONS AND 
PERSISTENCE OF HOMOCLINIC ORBITS 
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Abstract. We consider a dynamical system, possibly infinite dimensional or non- 
autonomous, with fast and slow time scales which is oscillatory with high frequencies 
in the fast directions. We first derive and justify the limit system of the slow variables. 
Assuming a steady state persists, we construct the stable, unstable, center- stable, 
center-unstable, and center manifolds of the steady state of a size of order 0(1) and 
give their leading order approximations. Finally, using these tools, we study the per- 
sistence of homoclinic solutions in this type of normally elliptic singular perturbation 
problems. 

1. Introduction 

A singular perturbation system usually involves different temporal or spatial scales. 
Here we focus on multiple time scales in which case the system takes the abstract form 
of 

(1.1) X = F{x, y, e) ey = G{x, y, e). 

The fast motions in the y direction are often some noise or transient behaviors and 
the slow motions in the x direction are more of the focus of the problem. In the 
singular limit as e — )■ 0, we obtain g{x, y, 0) = 0. Suppose y = 4){x) (without the loss 
of generality, assuming (/> = 0) solves this equation, the limit motion of x is given by 

(1.2) i; = F(x,0,,0). 
Let y = |, the y equation in (jl.ip takes the form 

/-I o\ ~ Gy{x, 0, 0) _ _ 

(1-3) y = — y + 9{x,y,e). 

€ 

The singular perturbation system (jl.ip is called normally hyperbolic if, for each x, 
the linear flow e**^^^^'^'*^^ on the y space is hyperbolic, i.e. it is exponentially con- 
tracting on one closed subspace and exponentially expanding in an complementary 
subspace. In this case, the standard normally hyperbolic invariant manifold theory 
\F1\ IHPSl IHej IBLZH IBLZ2] applies to yield a persistent normally hyperbolic invari- 
ant slow manifold given by a graph y = €(p{x,€). In the fast (and hyperbolic 
in natural) motions outside M^, solutions usually approach a neighborhood of Mg 
exponentially along its stable direction. After some time moving along the slow mani- 
fold, the solutions leave the neighborhood exponentially along the unstable directions. 
These motions of multiple scales can be connected by tools such as invariant foliations 
|F2t IF3t IHPSt ICLLt IBLZ3j and this geometric approach has led to a huge success 
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in the study of the dynamics of singular perturbation system p.ip . See for example 
[FillJ^ lJKllBLZi] . 

In the normally elliptic case, i.e. e*'^^^^''^''') is oscillatory instead of hyperbolic, on the 
one hand, the persistence of the slow manifold is not always guaranteed |GLH [0112] . 
On the other hand, solutions starting near {y = 0} should stay there at least for some 
0(1) time period due to the lack of strong exponential instability in the y direction. 
One typical situation of this type is when Gy{x, 0,0) is anti-self-adjoint. 

In this paper, with applications to both ODEs and PDEs in mind, we study these 
normally elliptic singular perturbation problems in an infinite dimensional dynamical 
system and possibly non- autonomous framework. Assuming Gy{x, 0, 0) = J, a constant 
anti-self-adjoint operator, we first justify the limit equation (|1.2p of the slow variable 
X through a careful averaging. 

A more important question is how much of the dynamical structure of the limit slow 
system (jl.2p remains in the singularly perturbed system (jl.ip . Elliptic type motions in 
the slow directions such as periodic or quasi-periodic solutions may be resonant with 
the oscillatory fast motions in the y direction. Some results have been obtained on the 
persistence of periodic orbits for nonresonant e < < 1 |GLH IGL21 ILol IMal ISZ2j . Here 
instead we focus on the basic hyperbolic structure - the local invariant manifolds near 
a steady state. Suppose (0, 0) persists as a steady state of (jl.ip for e << 1. Assume the 
linearization of the limit slow system (jl.2p has invariant stable, unstable, and center 
subspaces X'^'^''^. For the expanded system (jl.ip . the normal directions - the Y space 
- with the oscillatory linearized flow e*'' should obviously be considered as additional 
center directions. The first observation is, even though system (II. ip is singular, the 
existence of local invariant manifolds of (0, 0) is guaranteed by the standard theory (see, 
for example, |BJllHallCL] ) after a rescaling of the time by a factor of e. However, since 
the exponential growth/decay rates in the unstable/stable direction are 0(e) after the 
rescaling, this approach would only yield local invariant manifolds of the size of 0(e), 
which is far from being useful in most applications, such as studying the persistence 
of homoclinic orbits. 

Our main result in the manuscript is the existence and smoothness and the leading 
order approximation of invariant manifold of the steady state of the size of 0(1) based 
on a combination of the averaging and Lyapunov- Perron integral equation methods. 

As an application which is also an fundamental problem itself, suppose there exists 
a homoclinic orbits in the limit slow system (jl.2p and we study its persistence in the 
singular perturbation system (jl.ip which can be either weakly dissipative or conserva- 
tive. In the former, we derive the Melnikov function, which include an additional term 
coming from the fast directions, whose simple zero indicates a persistence homoclinic 
orbit to (0, 0). In the latter, when the system is analytic in reasonably low dimensions, 
along with some other structures such as the Hamiltonian setting or the reversibility, 
it has been shown that the stable and unstable manifold miss each other by an error 
like 0{e~~) |SunjlGetlLollTo] . Without these assumptions, we prove that there always 
exist orbits homoclinic to the center-manifold, forming a tube homoclinic to the center 
manifold. While we follow the well-developed geometric ideas in the finite or infinite 
dimensional regular perturbation problems |GHj IHMj ILMSWj ISZ3] , the proof heavily 
depends on the invariant manifolds we studied. 
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Before finishing the introduction, we would like to give two simple examples which 
partially motivated us to study this subject, while it is also easy to come up with 
examples in infinite dimensions. One is an elastic pendulum with fast and slow fre- 
quencies itself and the other one is a bifurcation problem which does not have any 
singular parameter in the appearance. 

A pendulum of the unit length with a fixed end is described by the Duffing equation. 
In a more careful model, the pendulum usually considered as rigid may have some small 
elasticity - meaning large elastic constant - allowing the pendulum to be stretched 
or contracted slightly in the radial direction. Let x be the angular and 1 + y be the 
radial coordinates, respectively, and the system takes of the form of a normally elliptic 
singular perturbation problem 



(1.4) 



(1 + y)x + 2xy + gsvux + 2e7(l + y)x - ^ ^ ^ Fi{x,y, e,t) = 
y - (1 + y)x'^ + \y - gcosx + 2ejy - eF2{x,y,e,t) = 



Where we also included the small damping and forcing. Formally, as e — ?• 0, i. e. the 
pendulum converges to be rigid, the corresponding singular limit ()1.2p for the above 
system becomes 

(1.5) y = 0, x + gsmx = 0. 

When there is no damping and the force is conservative, the problem is in the La- 
grangian setting and the limit equation is justified in [RU^ \Ai\ ITaj . In the dynamics, 
the state (vr, 0) is a hyperbolic steady state of (jl.Sp with a homoclinic orbit which of- 
ten leads to chaos even under small regular perturbation |GIIj . One may easily change 
the variables in the singular equation of y and make it anti-self-adjoint. Our general 
results apply to (jl.4p and give the criterion when the homolcinics persist under either 
dissipative or conservative perturbation. This example will be revisited in Section [6l 

The singular perturbations theory also applies to problems which may not be ex- 
plicitly in the form of (jl.ip . Consider an autonomous 4-dim ODE system with a 
parameter e which has the origin O as a fixed point for all e << 1. Assume, at e = 0, 
the linearized systems has simple eigenvalues iti and a double eigenvalue 0. While the 
unfolding of the focal point has been studied thoroughly (see, for example |CLWj ) . we 
note that the oscillatory motions are essentially at a much faster scale in the directions 
of the pair of elliptic eigenvalues. Under these assumptions, some simple normal forms 
transformations and near identity time rescaling, the generic form of the system looks 

::(:) ' m**) ^ ^ ''('^'^ ^ i"'^' * - ( -i ^l) ^ ^ «<i^i^ ^ i^i^' 

where x, y € and a«m(0) = 6(0) = 0. Rescale the system again by 

3 _i 
xi = exi, X2 = £2X2, y = ^y, t = e 

we obtain a singular perturbed systems in the form of (jl.ip of the normally elliptic 
type with the singular parameter = . If ^^^(0) > 0, the origin becomes hyperbolic 
in the x directions and we obtain the local center manifolds of order 0(1) size in the 
rescaled variables. If ^(0) 7^ in addition, we are in the right position to study the 
Hopf bifurcation from the eigenvalues itz in this rather degenerate case. (See |F5| for 
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an approach essentially different from the Hopf bifurcation.) A more detailed study of 
this type of bifurcation problems will be given in a forthcoming paper. 

The rest of the manuscript is organized as the following. In Section [T] we present 
the general framework and outline the main results on invariant manifolds and folia- 
tions. The justification of the limit slow equations and its linearization are obtained 
in Section [3j In Section [5] and [5] we study invariant manifolds and foliations and focus 
on their leading order approximations. Finally the homoclinic orbits are considered in 
Section [6l In the Appendix, we outline a process to block-diagonalize the linearized 
system of (jl.ip at a steady state. 

2. Framework and main results on invariant manifolds and foliations 

We formulate the problem as a non- autonomous infinite dimensional dynamical sys- 
tems with a singular parameter. This framework allows one to apply the general results 
to ODEs as well as PDEs or functional differential equations. 

In a Banach space Z and z G we usually denote a ball by Br{z, Z). Let Zi, Z2 
be Banach spaces and /c > 1 be an integer. We adopt the notations 

Lk{Zi, Z2) — L{(^^ Zi, Z2) = {bounded k — linear operators Zi ^ Z2] 

\(t>\- sup \(l){zi,...,Zk)\, for (p e Lk{Zi,Z2) 
ki|<i,-kfcl<i 

C^{Zi, Z2) = {h\h : Zi ^ Z2, fc-times countinuously differentiable with 
finite C'^ normj. 

Note Li{Zi, Z2) = L{Zi, Z2) is simply the space of bounded linear operators. 

Throughout the manuscript, we use D or to denote differentiations with respect 
to variables in the phase space and we will use d for derivatives with respect to time t 
or other parameters. 

Let X be a Banach space and Y a Hilbert space and we consider the system 




X = Ax + f{x,y,t,e) 
■ J 

y = -y + 9{x,y,t,e) 



The following assumptions may look complicated which is only due to our inten- 
tion to make the result applicable to PDEs where unbounded operators and different 
function spaces are involved. For ODE systems, these assumption would simply be 

• J is an anti-symmetric matrix and (/, g) are smooth functions. 

In general, we assume for some constants Co, 

(Al) A : Xi = D{A) — )• X, where Xi C X is endowed with the graph norm | ■ Ixu 
generates a Co-semigroup e*"^ on X such that je*"^! < Me'^^, t > 0, for some 
M > and w G M. 

(A2) J is an anti-self-adjoint operator on Y with domain D{J) = Yi, endowed with 
the graph norm \-\yi, which generates a unitary group e*"^. We further assume 

l"^" \l{Y,Yi) - ^0- 
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(A3) For A; > 1, 

(2.2) {D'f\D'g) G C^{Xi x Yi x M^L^(Xl x Yi,Xi x Yi)), 0<i<k, 

(2.3) {D'f,D'g) G C'^iXi x x 

L((X X Y) 0^-1 (Xi X Yi),X X Y)), l<i<k, 

whose norms are all bounded by Cq. 
(A4) I • |x € C'^(X\{0}, M"^), where M"^ denotes the set of positive numbers. 
(A5) dj G C^{Xi xYix M2,Xi), D^dJ G C^{Xi xYiX M^L(Xi,Xi)), D^dtf G 

C°(Xi X Yi X ]R^L(yi,X)), 5jc/ G C°{Xi xYiX M^F), DStg G C°{Xi xYiX 

R2,L(Xi X X y)) which are all bounded by Cq. 

Here the global boundedness are not important as we can always multiply them by a 
cut-off function. When unbounded, the linear operators often appear in the form of A, 

lA, or , etc. as coming from the linearization of PDEs [pj. The nonlinearities 

usually satisfy the above assumption on function spaces which are algebras or slightly 
better. 

Remark 2.1. In fact we can replace J by J(e) for each small e, then all results in this 
manuscript still hold except Proposition \4-i^ and \4.15 . 

Throughout the manuscript, C denotes a generic constant, possibly with subscripts, 
which may have different values as in different lines, and it only depends on the quan- 
tities involved in (Al)-(A5). Let C be another generic constant, possibly with sub- 
scripts, and the dependence will be specified in the context. 

When \y\ « 1, formally from ()1.2p . we expect the x equation can be approximated 
by the singular limit 

(2.4) xo = Axo + f{xo,0,t,0). 

However, in the normally elliptic singular perturbation problems, there is usually not 
a persistent slow manifold and we will first prove the convergence to (12. 4p . 

Almost invariant slow manifolds. In order to justify this limit, one need to es- 
timate the y equation in (12. ip . One of the key issues is to handle the 0(1) driving 
force g{x, 0, t, e) which occurs even at y = 0. It is very natural to first carry out a 
transformation 

(2.5) yi=y + eJ-^g{x,0,t,e) 
which yields 

(2.6) yi = ^yi+ gi{x,y,t,e) gi{x,y,t,e) = 0{e) + h{x,yi,t,€)yi. 

At y = 0, the driving force in this equation gi{x,0,t,€) = 0(e). One may repeat this 
procedure and obtain 

yk = ^yk + 0{e'') + 0{\yk\). 

Therefore this sequence of transformations yields an almost invariant slow manifold, 
close to {yk = 0}, with an error (to the invariance in the equations) of 0{e^). While this 
process increases the accuracy at the cost of the smoothness of the equation, it would 
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not give an invariant slow manifold and also, in our general setting, the unbounded 
operators A and J without other assumptions could bring other complications. We will 
work directly with p. II) in most part of the manuscript. However in this manuscript, 
with these transformations in mind, we actually often prove estimates with upper 
bounds in terms of g{x,0,t,€), so that they would yield much finer estimates when 
combined with a sequence of transformations as in the above. For example, see Lemma 
13.11 Lemma 13.31 Proposition 14.121 Remark 16.61 to see statements of this type. 

Theorem 2.2. Assume (Al) - (A3) and (A5). For any to e R, let T > and 
{x{t),y{t)), xo(t) be solutions of (j2.ip and (|2.4p on [to,to + T]. Suppose 

\x{to) - xo(to)|xi + \y{to)\Yi < Cie, 

there exists a constant C' which depends on M,uj,T,Co,Ci, and |xo(io)|xij such that 
then for any t € [to,to + T], 

(2.7) |x(t)-xo(t)|xi + <C'e. 

A more careful estimate of approximations can be found in Lemma |3. 11 In addition 
to the convergence of solutions on finite time interval, we also need the convergence of 
solutions of the linearized equations. Linearize ()2.ip and we obtain 

' 6x = A6x + Dxf{x, y, t, e)dx + Dyf{x, y, t, e)5y 



(2. 



5y = ^Sy + D^g{x, y, t, e)Sx + Dyg{x, y, t, e)5y, 



Let ^{t^tQ^x^y^e) be the solution map of (|2.ip . From the above equations assump- 
tions (Al) - (A3), and the Gronwall inequality, it is clear to see that is bounded 
uniformly in e. Higher order derivatives of <1> in x, y can be estimated in a similar 
way. In the leading order approximation of (12. Sp . we combine the linearized ()2.4p and 
a linearized y equation 

' (5x0 = A^xq + Dxf{xo, 0, t, 0)6x0 

(2-9) < . J 

SVo = -Syo + Dyg{xo, 0, t, 0)Syo. 
V e 

Theorem 2.3. Assume (Al) - (A3) and (A5) for k > 2. Let {5x{t),5y{t)) and 
{5xo(t),6yo(t)) be solutions of (|2.8p and (|2.9p . respectively. Suppose 

(2.10)x(io) - xo{to)\x, + \y{to)\Y, < Cie, 

(2.1lK|feo(to)|xi + \Syo{to)\Y,) + \6x{to) - 6xo(to)\x + \Sy{to) - 6yo{to)\Y, < Cie. 
Then there exists a constant C' which depends on M,u,T,Co,Ci, \x{to)\xi, such that 
\6x{t) - 6xo{t)\x + \Sy{t) - 6yo{t)\Y^ < C'e 

for all t G [to,to + T]. 

We can not obtain the estimates on \5x{t) — (5xo(t)|xi even if we assume \5x{to) — 
(5xo(to)|xi < Cie unless both |(5yo(io)|yi ^ Ci^- See Lemma [3.3^ Remark 13. 4^ and 
Remark 13.51 These theorems will be proved in Section [3l 

To study the local invariant manifolds, suppose (0, 0) is always a steady state and 
the limit systems is autonomous, i.e. 

dtf{x, y, t, 0) = dtg{x, y, t,0)=0 /(O, 0, t, e) = giO, 0, t, e) = 0. 
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We assume linearized (|2.4|) at has the exponential trichotomy, i.e. there exist closed 
subspaces X"'*''^ such that there exist constants oi < min{a2,0} and a'2 > max{0,a'^} 
and for t > 0, 

|ei(^+/.{0)) < C^^ait |g-t(A+/.(0)) < ^^g-a^t 

|ei(A+/.{0))|^^ < c^^a[t |g-t{A+/.(o))|^^ < Cie-"''. 

Moreover, we assume the linearized flow e*'-"^^^^'^^-' satisfies the same assumption as 
g-i(A+/a;(o)) and thus the expanded center space of (j2.ip should be X'^ © Y. Along 
with a few other technical assumptions, rough our main results on invariant manifolds 
and foliations in the phase space Xi x Yi is 

Theorem 2.4. For e << 1, in the space Xi x Yi, 

(1) There exists smooth invariant stable, unstable, center- stable, center-unstable, 
and center integral manifolds of (0, 0) which can be written as graphs of smooth 
mappings from a 6 -neighborhood of the corresponding subspaces to the comple- 
ments whose norms and 5 are independent of e. Moveover their derivatives 
in to, the time parameter of integral manifolds, is 0(e) when evaluated in the 
norm | • |x + | • |y- 

(2) The center-stable and center-unstable manifolds are foliated into the disjoint 
union of smooth families of smooth stable and unstable fibers which also written 
as graphs of mappings whose norms are bounded independent of e. 

(3) The stable and unstable manifolds are 0{e) close to those of (j2.4p . 

(4) The center- stable, center-unstable, and the center manifolds at {y = 0} are 
0(e) close to those of ()2.4p and their tangent spaces there are 0(e) close to the 
direct sum of the unperturbed ones and Y , respectively. 

Here by the term an integral manifold, we mean a family of manifold M(t) param- 
eterized by t so that the solution map of (12. ip starting at initial time to and ending 
at ti maps M(to) into M(ti). They are independent of t if the system is autonomous. 
The precise statement of these results of the invariant manifolds are given in Section 
Hand [51 

3. The singular limit system on finite time intervals 

The basic idea to handle the singular terms in the proofs of Theorem 12.21 and 12 . 3 1 is to 
average in time which appears in the estimate as integration by parts. Instead of ()2.4p . 
we consider the following regular perturbation problem as an initial approximation 

(3.1) x* = Ax* + /(x*, 0, t, e). 
In the rest of this section, we will use the notation 

(3.2) gl(x,t)=g(x,0,t,e). 

Lemma 3.1. Assume (Al) - (A3) and (A5). For any to eR, letT > and (x(t),y(t)) 
and x*(t) be solutions of (j2.ip and ()3.ip on [to, to + T] such that x(to) = x*(to) 
and y(tQ) = 0. Then there exists a constant C which depends on M,uj,T,Cq, and 
\x*(to)\xi, such that for any t € [tQ,tQ + T], 

\x(t) -x*(t)|xi + \y(t)\Yi < C'e(\gl\(joci(XixR,Y) + \Dx91\c"{XixR,l{x,y))) 
where C^C} denotes the space of functions in t and in x. 
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Proof. By (|2.1|) . (|2.4|) and variation of parameters formula 

(x - x,){t) = f e(*-^)^(/(x, y, r, e) - /(x., 0, r, e))dT, 
J to 

rt J 

(3.3) y{t) = e^^~'^^^{g{x,y,T,e) - g{x^,0,T,e) + g{x^,0,T,e)dT. 

Jto 



Due to the oscillatory nature of e* ^ , we integrate the last terms by parts 
/ e(*-^)Tg(x„0,r,e)dr = -e^'--^TeJ-'g{x.,0,T,e)\l 

+ / e(*-^)-eJ~^ {dtg{x^,0, r, e)dr + D^,g{xe,0, r, e) (Ax, + /(x,, 0, r, e)))(iT, 
-'to 

where we also use assumption (A3) to ensure the last term on the right hand side is 
well defined. Therefore, 
rt J 

/ e(*"^)"g(x*,0,r,e)dr < C e{\gl\cocl(Xixm,Y) + \Dx9l\co{XixR,L{x,Y))) , 

Jto ^1 

where C depends on Co and |x*|_Yi. Consequently, 

(3.4ij|x-x*|xi +|y|yi)(i) <C f (|x-x*|xi + \y\Y^){r)dT + C'e{\gl\coci(x^xR,Y) 



to 

+ \Dx9l lco(Xi xR,L{X,Y))) ■ 

Then the desired estimates follows from the Gronwall's inequality. □ 

Proof, of Theorem \ 2. 2\ Let (xi(i), yi(t)) be the solution of (j2.ip with the initial values 
2^1(^0) = 2;(to) and yi(to) = and x*(t) be the solution of p.ip such that x*(to) = x{tQ). 
On the one hand, from Lemma |3.H for any t € [to, to + 

\xi - x*|xi + \yi\Yi < C'e. 

On the other hand, by using the variation of parameter formula and the Gronwall's 
inequality, it is straight forward to show, for any t G [tQ,tQ + T], 

\x - xilxi + \y - yi\Yi + \x* - xolxi < C'e 
and thus the theorem follows. □ 

Remark 3.2. Combining Lemma \3.1\ with the iteration of the type of the transforma- 
tions ()2.5p . we may obtain asymptotic expansions of solutions of ()2.ip with the leading 
order term given by solutions of (13. ID and the error of 0{e^). 

For the linearization, we consider the following as the principle approximation 

5x^, = A5x^ + Dxf{x^,,0, t, e)5x* 
Sy* = i^y* + Dyg{x^ , 0, t, e)6y^ . 



(3.5) 



Lemma 3.3. Assume (Al) - (A3) and (A5) for k = 2 and use the same notations 
as in Lemma 13. it Let {5x{t),5y{t)) and ((^x*(t), (5y*(t)) be the solutions of (j2.8p and 
()3.5p respectively such that 

{5x{to),6y{to)) = {5x^{to),5y*{to)) and \5x(to)\xi + \Sy{to)\Yi < 1- 
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Then there exists a constant C depending on M,uj,T, and Co, such that for any 
t G [to,to + T], 



\Sy{t) - (5y*(t)|yi < C'ei^\gl\(.oci(^XixR,Y) + l-^^-f*lcicO(XixR,L(x,y)) 

+ \Dxdtgl\cO(XixR,L{Xi,Y))) ■ 

Plus, \6x{t) -6x^{t)\x^ < C'e i/ |(5y*(to)|yi < Cie. Otherwise \bx(t) - bx^{t)\x < C'e. 
Proof. By standard semigroup theory in Banach space and (A2), we have 

(3.6) (^\{6x,dy)\ + \{6x^,5y^)\^^ 

First we use (12.81) and (13.51) to obtain 



<C'. 

' CO{[to,to+T],XixYi) 



6x — (5x* = A{5x — 5x*) + [Dxf{x, y, t, e) - Dxf{x^,,0, t, e))6x^ 

+ Dxf{x, y, t, e){6x - 6x^) + Dyf{x, y, t, e){6y - 6y^) + Dyf{x, y, t, e)6y^ 

5y - Sy^ = (7 + Dyg{x^{t),0,t,e)){6y - 6y*) + Dxg{x*,0,t,e)5x 

+ {Dxg{x,y,t,e) - Dxg{x*,0,t,e))6x + {Dyg{x,y,t,e) - Dyg{x^,0,t, e))Sy. 
By using assumption (A3) and (j2.8p . we can write 
(3.7) 6x - = A{6x — 5x^) + hi{t, e) + Dyf{x, y, t, e)5y^, 

5y - 5y^ = (7 + Dy9{x*{t),0,t,e)){6y - 5y^) + /i2(t,e) + Dxg{x^,0,t,e)5x, 



where by (13. 6p . 

\hi\xi < Co{\Sx - Sx^lxj + \5y - (5y*|yJ + C'{\x - + blyj 
l^iU < Co{\5x - 5x^\x + \Sy - 5y*|y) + C"(|x - + blyj 
|/i2|yi < C'{\x-x^\x, + lylyj. 
In the rest of the proof, we will simply write 

Dygl (t) ^ Dygl {x, (t) , t) = Dyg{x, (t) ,0,t,e) 

and similarly for /, g, and gx etc. By assumptions (A2) and (A3), - + gy generates an 
evolution operator i?*(t, s) or £'^,(t, s; x*(s), e) which satisfies for t>s, 

dsE{t,s) = -E{t,s){:^ + Dygl{s)) |i?U(y,y) + \E\l(y,,y,) < Coe^"^'-'\ 

The most troublesome term in the second equation is gxSx. We use ()3.6p and integrate 
by parts to obtain 

(3.k)/ E{t,T)Dxgl{T)6x{T) dr = e E{t,to){J + eDygl{to))-^ Dxgl{to)6x{to) 
J to ^1 

-{J + eDygl{t))-^DxgmSx{t) + [ E{t,T)dMj + eDygl)-^Dxgl5x)dT) 

Jtn ^ / n 
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Prom assumptions (A3) and (A5) and Theorem 12.21 and equations (|2.8|) and (|3.6|) . we 
obtain 

ft 



[ E{t,T)D^gl{T)6x{T)dT < C'e(\D.^gl\co{XixRMX,Y)) 

which, along with the variation of parameter formula, the estimate on /i2, and the 
Gronwall inequality, implies the desired estimate on 5y — dy.^ . 

When |(5*(io)|yi ^ Ci^-, it is clear |5y*(t)|yi < C'e and thus the estimates on 5x — 5x^ 
also easily follows from the variation of parameter formula, the estimate on hi, and the 
Gronwall inequality. Otherwise, to deal with the most trouble term Dyf{x,y,t,e)5yif 
in the x equation, we use p.Sp to write 

6y,{t) = eJ~^6y,{t) - eJ^^ gy{t)5y,{t). 

Integrating by parts, we obtain from (A5) and ()3.6p 

I f e('--^^Dyf{x,y,T,e)5y.dT\ 



'to 

r-t 



C'e(l+\ dr{e^'-^^^Dyf{x,y,T,e))J~Hy,dT\^) < C'e 

and thus the estimate on 5x — 5x^, follows from the variation of parameter formula, the 
estimate on /ii, and the Gronwall inequality. □ 

Remark 3.4. If Sy^,{tQ) = 0, it is clear that \5x — 5x^\xi satisfies the estimate as 
5y — 6y.f. Without the assumption |5y*(to)|yi < C\e, we can not obtain the estimate 
on \5x — (5xo|xi since in the last step of integration by parts, there is a term 



f Ae^'-^^'^Dyfix, y, r, e)J-^6yodT, 

Jto 



which is only in X under current assumptions. 

Proof, of Theorem \2.3\ Let {5x^{t), 5y^{t)) be the solution of (j3.5p with initial value 
{5x{to),6y{tQ)). The estimates on (5x^, — Sxq and 6y^ — 6yo follow from the standard 
Gronwall inequality, which along with Lemma 13.31 implies Theorem 12. 3i □ 

Remark 3.5. Following from the same proof, if we assume instead (j2.1ip by 

\Sxoito)\x^ + \6yo{to)\Y, < 1 \Sx{to) - 6xo{to)\x + \Hto) - <5yo(io)|y < Cie. 

Then there exists C such that for all t € [tQ,tQ + T], 

\6x{t) - 6xo{t)\x + \6y{t) - 5yo(i)|y < C'e. 

4. Invariant Manifold 

In this section, we study the local integral manifold (as the system may be non- 
autonomous) of a stationary solution of (12. ip . namely, the center- unstable (stable) 
manifold, unstable (stable) manifold and etc in the framework of the Lyapunov-Perron 
integral equation. The main point is to obtain these manifolds of size 0(1) and their 
leading order approximations. Hypotheses (Al)-(A4) will be assumed and (A5) will 
be needed in some theorems with from Theorem 14.71 as specified. 
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4.1. Preliminary. In addition to (Al)-(A4) we assume 

(Bl) dtf{x,y,t,0) = dtg{x,y,t,0) = 0, 
(B2) f{0,0,t,e)=g{0,0,t,e) =0, 

(B3) When k = 1, assume {Df,Dg) are equicontinuous functions in x,y and e with 
respect to t at x = 0, y = 0, e = 0, i.e. for any s > 0, there exists S > such 
that if \x\x^ < 6, \y\Yi < |e| < for any t € R, 

\Df{x, y, t, e) - Df{0, 0, t, 0)\ l^x,xYMMXxY,x) < ^' 
\Dg{x,y,t,e) - Dg{0,0,t,0)\^^^^^y^ y^^^^^^^yy^ < s. 

We will write 

fx,y = D,,yf{0, 0, t, 0), g.,y ^ D^,yg{0, 0, t, 0) 

which are independent of t. 
For e Xi xYi, let 

(4.1) Fi{x,y,t,e) = f - fxx - fyy Gi{x,y,t,e) = g - gxX - gyy 
and A(r) € C^{R) such that 

A(r) = |^' , |A'(r)l<3. 

[O, r>l 

To take the advantage of the linear approximation, we cut off the nonlinearity 

(4.2) Fix,yXe) = A(^^^i±^)Fi G(x,y,t,6) = A(^^^^±^)Gi, 

r r 

then by assumption (A3), (A5) and (B2), F and G satisfy: 

F(0,0,t,e) = G(0,0,t,e) =0 \F\x, + \G{x,y,t,e)\Y, <rr, 

^ ■ ^ |Z)(F,G)|i(Xixyi) + l^(^'G)U(^xy) <r 

where r = r (r, e) = Gr with C depending only on Co and 

, , f=f(r,eo)= sup |L'(Fi,Gi)|i(Xixyi) + l-C>(Fi,Gi)|i(xxy)- 
(4.4) \x\x^+\y\Y^<r, 

ee[0,eo), teM 

Clearly, (A3) implies lim r = 0. Let 

r,eo— ^0 

Af = A + U. 

Since (12. ID is non- autonomous, in order to construct the local integral manifolds at 
time we translate the equation by to ^-nd modify it by cutting off the nonlinearity 

' x{t) = Afx + fyy + F{x, y, t + to, e) 



(4.5) 



yit) = (7 + 9y)y + 9xX + G{x, y, t + to, e). 



Clearly the system is unchanged in the |-neighborhood of (0, 0) and we will construct 
its global integral manifold which can be characterized by the exponential decay as 
t —7- ±00 at a rate close to that of e^^f Ix^ ^"- Naturally we need the following weighted 
continuous function spaces. In general, given r/ G M and Z a Banach space, let 

C^{Z) = {ze C°(M^,Z)| sup e-'>^\z{t)\z < +00} 

^ ±£>0 J 
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with norm | • |^^^ ^, where is a parameter 

'' ' ±i>o 

In order to hand the hnear terms fy and in system (|4.5|) which are not there in the 
usual study of local invariant manifold, we will introduce the following spaces which 
allow us to average. Let 

Bt;{p) = {{x,y) GC±(Xi) X C7±(yi 

where we also use B^{oo) to denote the corresponding Banach space. 

4.2. Existence and smoothness of integral manifolds. We start with the center- 
unstable and stable manifold. Assume the exponential dichotomy of the linearized 
system at (0, 0) 

(B4) There exists a pair of continuous projections {Pg^Pcu) on X, such that Pg + 
Pen = Ix 0, clearly X = PgX ® P.^X and 

(4.7) X^'^^ = Pg^cuX are positively invariant under e*^-^, 

e^^f can be extended to a group on X'^'^. 

(B5) There exist constants ai < and ai < a2, such that 

{e^'^fPsxlx < Ke"i*|x|x for t>0, x€ X, 

\e^^fPcux\x + \e^^^^^^h\Y < Ke^'^Wxlx + \y\Y) for t<0, xeX,yeY 

Remark 4.1. Let P,Xi = Xf and PcuXi = Xf, (B4) and (B5) imply e*^/ and 
gti-+9y) satisfy the same estimates with all norms replaced by \ • \xi and \-\yi- 

Theorem 4.2. Assume hypotheses (A1)-(A4), (B1)-(B5), and there exists rj such 
that r],krj € (01,02) then 

(1) There exist r > and cq > and a mapping kg : Br{0, Xf^ x Yi) x R x (0, eo) 
Br{^,Xl) such that the family of graphs 7W^"(to) form a locally invariant 
center-unstable integral manifold o/(0,0) of (|2.ip . 

(2) A backward flow is well-defined on the center-unstable manifold. 

(3) hs{(,cu,^,y,to,e) is in ^cu o,nd and continuous in to with the norms inde- 
pendent of e. 

(4) If we assume hypothesis (A5') to be introduced before Proposition \4. 7[ then 

dtohs{;;e) G C°(S,(0,Xr X Y,) x M, X') 
with the norm independent of e. 



Throughout this manuscript, we will use notations Ix, Iy for identity maps on X, Y , respectively. 
With slight abuse of notation, we also use the projections and P^u to denote their composition with 
the projection from X x Y to X. 
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Similar statements are given in Theorem 14.91 I4.1UI and 14.111 for stable, unstable, 
center-unstable, and center manifolds. These invariant manifolds are constructed in 
small, of 0(1) though, neighborhoods of 0. Sometimes, we do want to track the 
invariant manifolds in larger ranges. This can be achieved by combining the local 
invariant theorems and Theorem 12.21 and 12.31 See Proposition 14.151 

As in the standard approach, we will work on the global center-unstable manifold 
of (j4.5p which yields the local center- unstable manifold of (|2.ip . 

For any r] G (01,02), there exists e > such that for any € [0,e), there exist 
ro, eo > satisfying that for any e € [0, eo) and r € (0, tq), it holds 

(4.8) a{r]) = min{(Ti, (72(77), ^3(7?)} > 0, 
where 

(4.9) ai = l-C^e, ^2(7?) = 1 - 3(— ^ + — ^ + 1) (r + Coe^) , 

aM = 1 - -(-^ + K + l){r + 2Cle). 
02 - 

In the rest of this whole section, we always assume ()4.8p . 
Given to £ to simplify notation, we write 

F{x, y,T + to, e) = F{x{T),y{T),T + to, e) = Fix, y,T + to, e) fyy = fyyir) 

and such notation also applies to G and g. Let 

e*^/ 



and for any {x,y) G B~{p) and ^ = {£,cu,Cy) G ^1" ^ 

[A.mUx, y, to, em = U{t, e)i + f U{t - r, e) f (f"^"^' ^ + + ^J^^/) \ 

Jo \ G{x,y,T + to,e) + QxX J 

+ f C/(t-r,.)f^^(^(-'^'^ + *°'^) + ^^^)')rfr. 

It is standard to verify that {x, y) G B~ (p) is a fixed point of £^ if and only if {x{t),y{t) 
is a solution of (j4.5p with {I — Ps){x,y){0) = {S,cu,S,y)- Thus we focus on the fixed point 
equation 

(4.11) ix,y) = ,%u{x,y,C,to,e). 

Lemma 4.3. For any r] with oi < 7/ < 02 and e+,r, eo satisfy (j4.8p . there exists po 
depending on \^cu\xi,\^y\Yi, K,e-^,a (a defined in ()4.8p ). such that for any e G [0, eo) 
and p G [po^oo], ^cu defines a contraction mapping on B~{p) under the norm \ ■ j",,^. 

Proof. To keep the exposition clean, we will skip in most places the parameters S^cu, 
S,y, to, and e, which are fixed in this lemma. It is easy to obtain from the definition. 



(4.12) supe-"* {Ps + Pcu)^cu{x,y){t) 



X 



< K\U\x, + + ^){\DF\co+e.\fy\)\{x,y)\l^ 
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Since {x,y) € (^i) x (Yi) and Af is a closed operator, one can verify {Pg + 
Pcu)^cu{x,y) G C^((-oo,0),X) and 



(4.13) 



Therefore, 
(4.14) 



{Ps + Pcu)-^'%u{x,y){t) 

--Af{Ps + Pcu)^cu{x,y){t) + F{x{t),y{t),t + to,e) + fyvit). 



sup e 



(p, + Pc«):^=^c«(a;,y)(t) 

at 



X 



< K\U\x, + 



K 



K 



+ l){\DF\co+e^\fy\)\{x,y)^ 



^a2-r] rj-ai 
Again, by the definition of and integration by parts 

(4.15) {I -Ps-Pcu)%u{x,y){t) 

(4.16) =e'^-^^3y)iy^ \ e'-'-^^^i+^y^G{x,y,T + to,e)dT 

Jo 

e Jo £ 



{t-T){J+egy)f:^ 



Consequently, we obtain 
1 



(4.17) 



sup — e 



iI-Ps-Pcu)^cuix,y){t) 



< 



+ ( + K + 1) 



\DG\co +2e\J-^\\g^ 



\{x,y)' 



a2-r] 

Using (|4.8|) . clearly, there exists po > determined by ICylyj, \(,cu\xi, K, e^, and a such 
that for any p G {pQ, +oo], the above inequalities imply that maps B~{p) to B~{p). 
To prove it is a contraction, we can estimate in a similar fashion 



sup e 



<( + 

0,2 — 11 r] — ai 



{Ps + Pcu){^cu{xi,yi) - ^cuix2,y2)){t) 

K 



)(|DF|co + e*|/j/|)(|xi - a;2l^,i,xi + l^i ~ ^/sl^^.^^yj, 



(4.18) 



sup e 



(P, + Pc«)^(=^c«(xi,yi) - Sr^u{x2,y2)){t) 



X 



^ + — ^ + 1) (l^^bo + (l^^i - x2\-^i^x, + lyi - y^\le.,Y,)^ 



a2-r] r]-ai 



sup e 



(/-P. -Pc«)(=^c«(xi,?/l) - ^cu{x2,y2)){t) 



Yi 



K ^ ^ . \DG\co + 2e\J-'\\g^ 



\{x2 -xi,y2- yi)|^_, 



02 - ^ e* 
Therefore, defines a contraction mapping on B~{p) under the norm | • |~^^. □ 

Many other proofs in this manuscript will be very much in the fashion of that of the 
above lemma in the sense that integration by parts often provides an effective way to 
provide an extra e in the estimate. We will skip some details. 
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For any {£,,to,€) € X^^ x Yi x M^, let be the fixed point of ,%u, and 

(4.19) hs{^,to,e)=Psx{0)= r e-^'^fPs(F{x,y,r + to,e) + fyy)dT, 

J — oo 

From the standard argument based on the uniqueness of the contraction, one can show 
that hs and the center-unstable manifold Mf^{tQ), for any to £ 1^ are independent of 
7] G (01,02). Moreover, the flow map of (j4.5p starting on Aif^{tQ) is well-defined both 
forwardly and backwardly in t and, from ti to t2, it maps A4f^{tQ + ti) to A4f^{t() + t2)- 

Remark 4.4. Note since Xf^ xYi^ ToAif^, we cannot prove hg is bounded. See the 
Appendix. Also, as usual, A^^"(to) depends on the cut-off and thus is not unique for 



Smoothness in ^cu and ^y. Let z{S,){t) = {x{t),y{t)) denote the fixed point of =^(C)- 
From (|4.1ip . formally, 

(4.20) D^ziO = i^,^) = {I -D,^,u{z,to,e))-^U 

where 

+ U(t - T) (P-{DF{Z,T + + Jyi>)\ 

with the parameter e skipped. It is easy to verify that (/ — D^^cuiz, to, e)) ^ and thus 
the right side of (j4.20p are well-defined with respect to the same exponential rate ry. 
In order to estimate z{S, + — z{^) — Dgz(^)^', consider 

(4.22) (/ - D,^,u{z, to, e)) {z{C + O - ^(0 " D^^OO 

= + eO - z{0 - D,,%u{z, to, e){z{^ + O - ziO) - Ui'. 

In the above right side the linear terms fy and Qx, which do not vanish at (0, 0), actually 
disappear. Moreover, the evolution operator has a uniform bound though it depends on 
e. Therefore the exactly standard argument |CH ICLL] (where no differentiation in t is 
need which would product \) implies the smoothness ol z{-,tQ,e) : Xf** x Yi — )• B~, 
for any rj' < rj. This establish the smoothness of Aif^{tQ) 

ii^ ^cu and ^y. Similarly, 

the smoothness of M'^ can be obtained if we assume that there exists 77 such that 
oi < r], kr] < 02. 

Dependence on to. In order to smoothness of A4f^ in to, in addition to hypotheses 
(A1)-(A4) for k = 1 and (B1)-(B3), we assume and 

(A5') {dtf, dtg) e CO(Xi xYiKm?,XxY), {Ddtf, Ddtg) G CO(Xi x ^ x M^, L{Xi x 
yi,A:xy)), {DdAf^DdAg) G C^{XixYixm?,L{XixYi,XxY)) Moreover, 
their norms are bounded by Co. 

Remark 4.5. In fact, the assumptions in (A5' ) on Ddt{f,g) and Ddedt{f,g) are only 
needed when one has to work with r] > 0. In fact, if r] < 0, our next theorem still holds 
if we only assume dt{f,g) G X x y dtde{f,g) ^ X xY. 
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To avoid dealing with x{t) as in the above contraction argument, which may intro- 
duce a factor of -, we introduce a shght variation of which will be somewhat more 
easily used in the proofs of the following proposition as no time derivative is directly 
involved. Given parameters ^ G Xf^ xYi, to, e, for z = {x,y) € C~{Xi) x C~(Yi), let 

(4.23) torn = Uit, e)^ + (0(t), ^(t)) 
where 

m = f e^'-^'^^f P,u{Fiz,T + to,e) + fyy)dT 
Jo 

+ [ e^'-^^^fPs{F{z,T + to,e) + fyy)dT 

J — oo 

and 

m = f e('--^(i+3y)G(z, T + to, e)dT - e( J + e^,)"^ {gMt) + e'^i+^y) g,x{0)) 
Jo 

+ e / [J + egy)-^e^'-^^^-+3y)g^{AfX + F{z,T + to,e) + fyy)dT. 
Jo 

We equip C~ {Xi ) x C~ {Yi ) with the norm 

(4.24) \\z\\\ = supe-''*(|x(t)U, + M^^), 

t<o (-* 



and we will also use || • ||.;; to denote 



-.'(l,(,)l, + Mf)k)^ 



= sup e 

t<o 

We denote the balls in these norms by 

(4.25) B~,{^) = [z\\\z\\\, S-(oo) = { 



z\ \\z\\^ < oo 



Obviously, is come up with from ^ after integrating by parts in the y component. 
Using (A3) , it is straight forward to prove that £^cu is still a contraction on {Xi ) x 
C~{Yi) under the norm in (I4.24p . Moreover, its linearization is also a contraction 
under both of the norms || • ||^ and || • ||^. Namely, for some < a' < 1 with a similar 
form as a defined in (|4.8p . we have 

(4.26) ^cu{z,to)--fcu{z',to) <{l-a')\\z-z'\[n \\D^uz\\rj < {1 - ct')\\z\\^. 

V 

By the uniqueness, and have the same fixed point. 

Remark 4.6. In the following, we will repeatedly use the fact that z{^), the fixed point 
of '%u Oind 3'cu, belongs to B~{oo) for any r] G (oi, 02) as long as (j4.8p is satisfied. 

Since the time derivatives involves unbounded operators A and J, we do not expect 
A1^m(*o) to be smooth in to in Xi xYi. 
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Proposition 4.7. Assume the conditions in Theorem \4-^ for k = 1, then 

/i,(-,-,e)GCO(Xr xyixM,Xf). 
// we further assume (A5' ), 

Remark 4.8. If in (A5) and (A5') we assume the smoothness of {f,g) : Xi xYi xM^ 
Xi X Yi, the same proof implies that dt^hg € Xf is continuous. See also Remark \4-17\ 

Proof. To simplify notations, we will also ignore the e variable in F, G. We first claim 
(4.27) lim ^cu{z,ti) - Wcu{z,to) ^ = 0, 

for any ai < rj < rj' < 02 and z € B~,{oo), where B~,{oo) = ^z\\\z\\^, < oo|. In 
fact, for any s > 0, we will show the above quantity is bounded by some Cs when 
\ti - to\ « 1. Let T2 = By (A3), DF{pz{t),t'), DG{pz{t),t') are uniformly 

continuous on {p,t,t') G [0, 1] x [T2,0] x [T2 + tQ — l,to + 1]. Therefore, there exists 
s' > such that if \ti — to| < s', 

\DiF, G){pz{t),t + ti) - D{F, G){pz{t), t + io)L(.Y,xn,Xixyo < « 
for {p,t) G [0,1] X [r2,0]. Rewrite F(z(t),t + ii) -F(z(t),t + to) as 
F{z{t),t + ti)- F{z{t),t + to) 
=F{<t).t + ti) - F{0, t + ti)- F{zit),t + to) + F{0, t + to) 

= (^j^ DF{pz{t),t + ti) - DF{pz{t),t + to)dp^ z{t). 

It follows that for |ti — to I < s' , 

(4.29) |F(z(t),t + ti) -F(z(t),t + to)|^^ < s\z{t)\x,xYr, 

which is also true for G. To obtain (j4.27p . we split the integration intervals in the 
definition of into t < T2 and t > T2. On the former, the estimate can be obtained 
by using (I4.29P and the exponential bound of z and on the latter we only need to notice 

\{F,G){z{t))\x, < Cr\z{T)\x,xY, < Cre^^ < Cfse^'\ for r > T2. 
An immediate consequence of (I4.27P is 

(4.30) ^„GC7°(5-(oo)xE,fl-(oo)). 

In fact, by the same procedure we can also prove the following stronger statements, 

(4.31) d¥,u G C\B~,{^) X M,L(5-,(oo),5-(oo)), 

where ai < rj < r]' < 02 and D is the differentiation with respect to z. 

Let Zi be the fixed point of 3^cu{-,ti) for i = 0,1. From ()4.26p it is easy to see 

\ZI - ZqWI^ < (1 - 0-')lkl - ZoWli + '%u{zo,ti) - ^cu{zo,to) 
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Together with (|4.27|) and Remark 14. 6t it imphes 

1 



(4.32) lira \\z{ti) - z{to)\\l, < ^ lim ^Uzo^h) - ,%uizo,to) 



= 0, 



From the definition of /is, we obtain its continuity in to- 

To prove the second part, by our assumptions in (A4) and (A5') involving X and 
y, one may also prove 

(4.33) d¥cu G C\B-,{^) X M,L(:B^;;,(oo),B,7(oo))) 

(4.34) dt,^,u G C0(5-(oo) X M,;b7 ((»)), 

where (oo) = < oo|. 

Since (B2) implies F(0, t) = G(0, t) = 0, we write 

F{zo,t + ti)- F{zQ,t + to)-dt,F{zQ,t + to){ti-to) 

(4.35) / Ddt,F{qzQ,t + pti + {l-p)to) 
- DdtF{qzo,t + to)dqdp) {zo{t)). 



Assumptions (A5'), (|4.35p . and a similar estimate for G yield \\dtQ^cu{z)\[r] < C"||z||^, 
which implies, for any z G B~,{oo), 

(4.36) ^cu{z,ti) - ^cu{z,to) <C'\ti-to\\\z\\l, 
where C depends on C,ai,a2,T]. Since 

(4.37) zi- zo = ^cu{zi,to) - ^cu{zo,to) + '%u{zi,ti) - ^cu{zi,to), 
using ()4.33p and ()4.36p . we obtain 

C 

(4.38) Iki -^ollr, < -^\\zi\\li\ti - to\. 
We continue to write zi — zq as 

(4.39) zi-zo = dto^cu{zo,to){ti - to) + D^cu(.zo, to){zi - zo) + Ri + R2, 
where 

(4.40) -Ri = ¥cu{zi,ti) - ¥cu{zi,to) - dto¥cu{zo,to)iti - to), 
R2 = '%u{zi,to) - ^cu{zo,to) - D.%u{zo,to){zi - Zo). 

By ()i:32D and gSlD, 

ll-Rillr? =\ti - to\\\ / d^3^cu{zi,pti + (1 -p)to) - dto'9'cu{zo,to)dp\\ = o{\ti - to\). 
Jo ' 

Using (]4.33p . we have |!i?2||»7 = o{\\zi — zo\\rj) = o{\ti — to\)- Therefore, by (I4.39p . 

(4.41) dt,zo= (^I - D^,u{zo,to)) ^ dt,-fcu{zo,to). 

and we obtain dt^hsi-, •, e) G C°(Xi^" x x R, X'). □ 
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Since the cut-off function does not change the system in a neighborhood of radius 
I and r is taken independent of e, we obtain Theorem 14.21 

From the exponential dichotomy, we can also construct the stable integral manifold. 
For G Xf and {x,y) G B^{oo), define 

6, to, em = Uit, e% + U{t - r, .) (^^ ^' ^ + + ^^^)) dr 

(4.42) +r C/(t-r,e)f^^-(^i^^'^'l+*°\l+^^^)V^' 

where Fi , Gi are introduced in ()4.ip . One may note that Fi and Gi are not cut off 
in ^ as opposed in the construction of the center-unstable integral manifold. This 
rather standard practice is due to the fact ai < and thus the local information near 
the steady state along is sufficient to determine the unique stable integral manifolds 
where the solutions decay exponentially. Following the same proof as in Lemma 14.31 
one can prove that if (|4.8p is satisfied, ,% defines a contraction on B^{p), where p is 
sufficiently small but independent of e << 1, under the norm | • |+^^ given in (14. 6p . 
Therefore 

Theorem 4.9. Assume hypotheses (Al)-(A^), (B1)-(B5), and there exists r] such 
that r],kr] S (01,02) then 

(1) There exist r > and eo > and a unique mapping hcu '■ Br{0,Xl) x R x 
(0, eo) Br{0, Xf" X Yi) such that the family of graphs TWf (to) form the locally 
invariant stable integral manifold of (0, 0) of (12. ip . 

(2) Solutions are on if and only if they decay with exponential rate r] as t ^ 
+00. 

(3) hcu{^s,tQ,e) is in and continuous in to with the norms independent of e. 

(4) If we assume hypothesis (A5'), then 

dtMr,e) € CO(i?,(0,Xf) x M, X™ x Y) 

with the norm independent of e. 

In the following, we will give the hypotheses on center-stable and unstable integral 
manifold. 

(CI) There exists a pair of continuous projections {Pcs,Pu) on X, such that Pes + 

Pu = Ix and X'^^''^ = Pcs^uX are positively invariant under e^^f . 
(C2) There exist constants 02 > 0, and a'l < a'2, 

le^'^fPcsxlx + \e^^^'^^^^y\Y < Ke<\\x\x + \y\Y) for t > 0, xeX,yeY 
\e^'^'Pux\x < Ke^'^^lxlx for t < 0, x£ X. 

As in Remark 14.11 e^^f and e^^~~^^^^ satisfy the same estimates with norms replaced 
by I • \xx and | • [vi- I^et Pcs^i = and PuXi = X". By the same proof. 

Theorem 4.10. Assume (Al)-(A4), (B1)~(B3), (C1)~(C2), and there exists rj such 
that r],kr] S (a'1,02) then 

(1) There exist r > and eo > and mappings 

hcs : Br{0,X^) X M X (0,eo) ^ Sr(0,Xf x Yi) 
hu : Br{0, X Yi) X M X (0, eo) ^ 5^(0, Xf) 
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such that the two families of graphs A4"(to) CLnd Al^*(to) form locally invariant 
unstable and center-stable integral manifolds o/(0, 0) of (|2.1|) . 

(2) A backward flow is well-defined on M."^ and solutions are on JvV^ if and only if 
they decay with exponential rate rj as t — oo. 

(3) hcs{^u,to, e) is in S,u and hu{^s,(,y,to,^) is in and S,y and both are 
continuous in tQ with the norms independent of e. 

(4) If we assume hypothesis (A5'), then 

dt,h,s{;e) G C°(5,(0,Xn X M, x Y) dt,hu{;e) G C''{Br{0,Xt x Yi) x M, X^) 
with the norms independent of e. 

Let Xf = Xf" n Xf = {Ix - Ps- Pu)Xi = PcXi. By taking the intersection of 
Aif^ and we can obtain a center manifold in the standard way. 

Theorem 4.11. Assume (Al)-(A4), (B1)-(B5), (C1)-(C2), and there exists rj± such 
that i]^,krj^ G (a'^^jflg) and r]-,kr]- G (01,02) then 

(1) There exist r > and eo > and mappings 

^ = {^s, -^u) : 5r(0, XI X n) X M X (0, eo) ^ 5^(0, x Xf) 

such that the family of graphs A^^(to) form a locally invariant center integral 
manifolds 0/ (0, 0) of (IZT]l . 

(2) A backward flow is well-defined on Ai^. 

(3) ^(^c; ^0, e) is in and and is continuous in tQ with the norms inde- 
pendent of €. 

(4) If we assume hypothesis (A5'), then 

9to^(-,e) G C°(5,.(0,X^ X Yi) X M, X* X X^ 
with the norms independent of e. 

4.3. Asymptotic estimates of Invariant Manifolds. In Theorems 12.21 and 12. 3^ we 

have demonstrated that p.4p can be viewed as the singular limit of ()2.ip as e — ?> 0. 
Therefore, one may expect the perturbed invariant manifolds should be close to the 
unperturbed ones. In this subsection, we will give the leading order approximation of 
the invariant manifolds. 

In this subsection, we will use the notation 

(4.43) Fl{x,t) = F{x,Q,t,€) Gl{x,t) = G{x,0,t,e). 

Instead of ()2.ip directly, we first consider compare (|4.5p with the following regular 
perturbation problem 

(4.44) X* = AfX^ F^{x^,t-^to) 

where we also included the dependence of F on e. Under assumptions in Theorem 14.21 
= is an steady solution and it has local integral manifolds. In fact, for ^cu £ X™ 
and X G C~(Xi), let 

X{x){t) =e'^fU+ f e^'~^'>^fPcuF:{x{T),T + to)dT 

(4.45) , ^° 

+ / e^'-^^'^fPsF:ix{T),T + to)dT. 

J —00 
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By the exponential dichotomy (B5), ■% is a contraction on {Xi). For S^cu £ ^T^^ 
let x^:{t) be the fixed point of Define 

(4.46) K{^,u,to,e) = PsX,{0)= / e-^^f PsF:{x.{t),t + to)dT 

J — oo 

MTito) = {U + KiCcu,to,e)\Ccu G 
which is the center- unstable manifold of (|4.44p . 

Proposition 4.12. Assume the hypotheses in Theorem \4-^ (^5) and (A5'), then 
\hsiCcu,0,to,e) - h*{£,cu,to,e)\x^ ^ C'e(^\gx\L{x,Y) 

+ \DxGl\c(\XixR,L{X,Y)) + l^t^*lcoci(XixR,y)) 
where kg is defined in (|4.19p . If k >2, we also have 
\^Cy^s{^cu,0,to,e)\^^y^^^,^^ < C'e 

\D(^^hs{£,cu,0,to,e) - D^^^h*g{^cu-,to, e) \ ^^^cu ^xi) 

<C'^{\9x\LiX,Y) + \^xGl\cOci(XixR,L{X,Y)) + \^t^l\c°Cl{XixR,Y)) 

where C depends on Cq, K,ai,a2, and \^cu\xi o,nd the norm C^C^ means in t and 
in X € Xi. 

The reason the estimate on D^^hg is only in Y is similar to that for Lemma 13.31 

Remark 4.13. The exactly same estimates in this subsection also hold for the center, 
center- stable, stable, and unstable manifolds except for the latter two, there is no D^^ 
involved. 

Remark 4.14. The reason we include some these terms in the above upper bound 
goes back to transformation (j2.5p . Conceptually, with certain (mild) assumptions on 
A and J, one may carry out a sequence of transformations in the form of (j2.5p to 
make gx = and G(x,0,t, 0) = 0(e'^^^) when measured in appropriate norms. There- 
fore the above estimates immediately implies that the integral manifolds of ()4.44p are 
approximations of those of (|4.5p at {y = 0} with an error of 0{e^). Since (|4.44p is a 
regular perturbation problems of (|2.4p . one may compute the Taylor expansions of the 
integral manifolds of (j4.44p up to the order 0{e^~^) and thus they also serve as the 
leading order expansions of the integral manifolds of ()4.5p at {y = 0}. 

Proof. We will denote Px and Py the projection to X and Y. Let (x, y) be the fixed 
point of c%u with parameters ^cu and = and x* be the fixed point of with the 
parameter ^cu- In the rest of the proof, we will use ,%u{x,y,e) to denote S^cu{z,to), 
which is introduced in (j4.23p . From (|4.26p . 

(4.47) \\{x-x,,y)\\\ < \\Ku{x,y) - ^„(x„ 0)||^ + ||^„(x*, 0) - ^(x,)||^ 

< (1 -a')||(x-x*,y)|l^ + ||=^„(x,,0) - =^(xo)|l^, 

which implies 

\\{x - x.,y)\\^^ < ^||^n(x.,0) - X{x,)\\l,. 
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From the definitions of S^cu and one may compute by integrating by parts 

(4.48) (^«(rE„0) - %{x,)){t) = Py^„(x*,0)(t) 

(4.49) = e(J + egyY^ ( e^^--)^-.^^y) (dtG% + (5. + I^xGj) + F,^))dr 

+e*(7+9«) (^^x* (0) + G% (x, (0), to)) - (ffx^;* (t) + (x, (t), to + t)) 



where 9iG^, DxG%, and F^*^ in the above integral are evaluated at (x*(r),to + T). Using 
assumption (B2), we immediately obtain 



(4.50) \\{x -x^,y)\\]^ <C'e{\gx\L{xx) + !^:r-GJ|c"'(XixR,L{x,y)) + \dtGl\coci(x^ 



xR,y); 



where we need \dtG%\(jO(ji to bound dtGl by and |x| provides the necessary 

decay in t. Consequently the estimate on hs — h* follows. 

To prove the second part, choose r/ such that oi < r],2r] < S2- Let {4>'^ (t) , tp'^ (t)) be 
the derivative of (x(t),y(t)) with respect to at {Ccui^) and (p*{t) be the derivative 
of x^,(t) with respect to S.cu, so we have 



(4.51) icj)\r) = D.%uix,y)icP',r) + e'^f , = D5;(x, )(</.*) + e*^/. 

As in Theorem 14.21 It is easy to show ||j;^(-js(^cii (Xi)) bounded uniformly in e for 
i = 1,2. By gSH), 



(4.52) 



L{X1",B-^{^)) 



+ 



L(XJ",B2- (oo)) 

(Z)^„(x,y) -Z?^„(x„O))(0*,O) 



L(Xf",B2-{oo)) 
L(X™,B2-{oo))' 



In the first term on the right side, D^cu{x, y) is bounded by 1 — a' according to (j4.26p . 
Using ()4.50p and the fact that F and G are C^, we obtain through straight forward 
computation 



(4.53) 



L(Xf-,B2-(oo)) 

<C' e[\gx\L(X,Y) + \DxGl\cO(XixW.,L(X,Y)) + l'9tG'^lcOci(XixR,Y)) ■ 

For the last term in ()4.52p . one may calculate 

DWUx*M<t>\^) - DX{x,)4>* = PYDWcu{x,M<t>\^) 
=e{J + egy)-^[ J^e^'-^'^^i^Sy) ^D^dtGl + DIGUAjx, + F^) 
+ {gx + DxGDiAf + DxFl))(t>,dT 

{gx + DxGl{x,{0),to))(l),{0) - {gx + D^G^(x,(t), to + t))Mif 



+ e 
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where D.dtGl, D^Gl, and D xF^ in the above integral are evaluated at (x=f(T),to + 
r). Along with <g^, (jt^ . and (gSSD, it implies 



+ |-C'xG^lcf'Ci(XixR,L(X,y)) + \^tGl\cOci{XixR,Y)) 



and thus the estimates on D^^^hg — D^^^h*. 

Finally, with slight abuse of notation, we still use (t/)^,^'^) to denote the derivative 
of with respect to at = 0. Using (|4.26p and (j4.50p . it is straight forward to 

show that, at ^y = 0, 



(4-54) \y\c-^y^^ < C'e |(0^ ^^)li(y,,5-{oo)) < C 

where r]' can be taken in a compact subinterval of (oi, 02). Clearly 

Like one can show PxDx^%uix,y) is a linear contraction on C^{X). Therefore 

(4.55) l'^lL(yi,c-,(x)) < C'\PxDyXix,y)r\LiY„c,'^{X))- 

To estimate the right side, we notice that satisfies 

i^tit) = {^ + 9y + DyG{x, y, t + to, e))ip%t) + {DxG{x, y, t + to, e) + 5x)<^'(t) 
which can be rewritten as 

r = eJ-^^l - eJ-\{gy + DyG)^ + {D^G + g^W). 
Substitute this identity into ()4.55p and use ()4.54p . we obtain 

\^'\l(YuC^^{X)) < C'e{l + \PxDy,%ix,y)J-^^Pt\L{YuC^^{X))) 
<C"e(l + supe-2^* ( [' Pcu+ r Ps^e^'-^^'^f {DyF{x,y,T + to,e) + fy) 

^ t<0 Jo J~oo ' 

J-^ll^fdT ] 

L(yi,x=")/ 

Integrating by parts and using ()4.54p to control y, we obtain the desired estimates. 
The estimates can not be improved to the norm Xf as A is produced in the integration 
by parts. □ 

To consider invariant manifolds in larger ranges, let $(T, to)-2,e) and ^*{T,to,x,€) 
be solutions of (j4.5p and (j4.44p from time to T — to (so from time to to T for (j2.ip 
and (|3.ip ) with ^{to,to, z,e) = z = {x,y) and ^*{to,to,x,e)x. We will skip writing 
to and e in <I> and <&* in the next proposition as it would not be altered. Combining 
Lemma [3. H Lemma [3. 3 1 Remark l3.4l 13.51 and Proposition 14. 121 we obtain in a straight 
forward manner 
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Proposition 4.15. If the hypotheses in Theorem \4.^ (^5) ^"^d (A5') hold for k = 2, 
then there exists C which depends on C, if, ry, ai, 02, r, e^, \T — t^], |^c«|xi such that 



XixYi 



(4.56) HT,C,u + hs{U,0)) - <^*(.T,U + KiU)) 
(4.5'^ C'e(^\gx\i(^x,Y) + \DxGI\co(XixR,l(x,y)) + l'9tG'J|co(XixiR,Y) j; 

(4.58) D^^AHT,Uu + h,{^,u,0)))-D^^J<^>*iT,U, + h:iC,u))) ^ 

L(X-^ ,XiXYi) 

(4.59) +\Py (% U + ^y + hsiU, ^y))) - E{T, U + KiU))) 

(4.6(g C'e{\gx\L{X,Y) + l-^a;G^lcOCi(XixIR,L{X,y)) + \9tGl\c0c^(^XixR,Y) 



Cy=0 

) 



L{Yi,Yi) 



LiYi,X) 



<G'e 



where E{T, to,^cu+hl{S,cu), e) is the evolution operator generated by j+Dygl{^*{t, ^cu+ 
^si^cu)),t) with initial time to and terminal time T and Px,Py denote the projection 
from X X Y to X and Y , respectively. 

See Remark 14. 141 for the explanation why the above upper bounds are taken in such 
a tedious form. Also, the reason PxD^^^ is estimated only in X is Lemma l3.3[ 

Since the system is autonomous when e = 0, we also expect the derivatives of the 
integral manifolds in to is of order 0(e). This will be used in studying homoclinic 
orbits. 



Proposition 4.16. Assume the same condition as in Proposition \4- 7| for k = 2 and 

iCyWi < Cie, then 

\dtohs{Ccu,(,y, •,(^)\co{K,x^) < c'e, 
where C' depends on Ci, |^cu|xi CLnd constants in assumptions. 

Remark 4.17. // in (A5) and (A5') we assume the smoothness of {f,g) : Xi x Yi x 
—7- Xi X Yi, the same proof implies dtahg € X^ is of 0{e). See also Remark \4.8[ 

Proof Let zo = (xo,yo) be the fixed point^of .%ui-,to) and (0,-0) = idtoXo,dtoyo). 
Prom (j4.4ip . we only need to estimate dtf)^cuizo,to). Notice, from assumptions (Bl) 
and (B2), 

Ft{z,t,e) = e Fteiz,t,Tie)dTi = e / DFte{T2Z,t,Ti)zdT2dTi 
Jo Jo Jo 

and similar estimate holds for G. These estimates immediately implies 

\\dto-!^cu{zO,to)\\2r] < C''e|[zo||2r7 

and the thus the proposition follows. 



□ 



5. Invariant Foliation 



With the center-stable (center-unstable) integral manifolds constructed in the 
previous section, we will give the sketch of the construction of the stable (unstable) 
fibres inside the center-stable (center-unstable) manifold under the same assumptions 
in this section. We will use the stable fibres as an illustration and similar construction 
also works for unstable fibers. 



NORMALLY ELLIPTIC SINGULAR PERTURBATIONS 25 

For icy = iiciy) G XI X Yi, let {x{icy){t),y{icy){t)) be the solution of (HSJ with 
the initial value (at t = 0) on the center manifold 

(5.1) i = icy + ^siicy: to, e) + ^u{Ly,to, e). 

The solution stays on ^A'^{t) and satisfies 

(5.2) = U{t, e)C + r U{t -r,e)( If^ ^' " + % | dr. 
\yit)j Jo ' \G{x,y,T + to,e) + gccxj 

To simplify our notation, for {x,y) € Xi x Yi, we write 

(5.3) F{x, y, icy, t, e) = F{x{icy){t) + x, y{icy){t) +y,t + to, e) 

-Fix{icy){t),y{icy){t),t + to,e) 

or very often in short as F{x,y,icy)- Such notation also applies to G. 

For each triple {is,ic,iy) & Xf x x Yi and ai < r] < 02, it is the standard 
knowledge that {x{t),y{t)) is a solution of (j4.5p satisfying 

P,(x(0) - x(^ej(0)) = 6 and (i,y) ^ (x, y) - ((x(^,j,), y(ecj;)) G 5+(oo) 

where B^{oo) was defined in (|4.6|) . if and only if (5?(.),y(.)) is a fixed point of 

J +00 \ '~^{X,y,t,cy) -r dxX J 

One first notices that, for fixed icy, has the same form as with only an 
additional parameter icy Moreover, by (|4.10p 

F{Q,Q,icy,e) = G{Q,Q,icy,e) = Q \DF\c- = \DF\c<^ <r \DG\c^. = \DG\co <r 

where D is the differentiation with respect to (5;, y) or (rr, y). Through exactly the same 
procedure as in Section U we obtain that 'Sg defines a contraction on B^{oo) under 
the norm | . 1+^^ defined in (|4.6p . Clearly, if is = 0, (x, y) = (0, 0) is the unique fixed 
point of (|5.4p . Moreover, in the study of the G^ smoothness of the fixed point {x,y) 
with respect to is, the linear terms fy and Qx, which are not small, disappear again. 
The evolution operator U has a uniform bound though it depends on e. Therefore the 
exactly standard argument in regular perturbations |CLLj (where no differentiation in t 
is need which would product i) applies and yield the smoothness of {x, y) € i?^(oo) in 
is. Therefore, by fixing e+ small and then choosing r and e sufficiently small accordingly, 
we obtain the following theorem 

Theorem 5.1. Assume (Al)-(A4), (Bl) - (B5), (CI), and (C2). If there exists rj < 
with ai < krj < rj < 02, then for each triple {is,ic,iy) G -'^i x Xf x Yi, (15. 4p has a 
unique fixed point {x, y) € -B^ (00) such that 

i) lfis = 0, ix,y) = iO,0). 

ii) {Dlx, Dly) G CO [xi x x Yi, 5+^(oo)) , where j = !,■■■ ,k. 
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Let {x,y) be the fixed point of (|5.4|) corresponding to and ^cy, we define 

(5.5) acu{^s,^cy,to) = e + (S(0),y(0)) W!{^cy,to) = {acu{Cs,Ccy,to)\Cs e ^l'} 

where is given in ()5.ip . Usually V^^ i^cy , to) is called the stable fiber based at ^. 

Remark 5.2. Clearly, {x + x,y + y){0) G A^^''(io), where {x,y) is the solution of (|5.2p 
with parameters (CcCy^'to) <^f^d thus W|(^cyi*o) C A^^''(t) 

To study the smooth of (^cj/j^o) with respect to £,cy, for a positive integer k >2, 
define 

Afe = G M^|ai < /c?? < ?7 < min{0, 02}, a'^ < j-q' < a^, 

(5.6) ai < T] + jr]' < 02, i = 1,2, ••• 

Theorem 5.3. For k > 2, assume (A1)-(A4), (B1)-(B5), (C1)-(C2), and Afc is 
nonempty. For any compact subset S of K^, by fixing e^, small and then choosing r 
and e sufficiently small accordingly, then for any {r],r]') € S, (|5.4p has a unique fixed 
point {x,y) G B^{oo) such that 

i) {Dij,Dij) G (Xf X XI X yi,i?+_^^,^,(oo)) , 

ii) (D'T^^dI X, DT'-^Di y)eC^( Xf x xYi,B+ _^ . , (00)) , 
where m = 2, • • • , /c, j = 1, • • • , m — 1. 

Our spectral gap assumption on is essentially the same as in |CLLj and the proof 
of the theorem again follows from the same procedure which is based on the definition 
of the Frechet derivatives. One may notice that Ocu was only proved to belong to 
'^h^ while we have it in ^^^cy if J < ^- In fact, it is easy to verify that the 
same proof works to yields our above slightly stronger version. For details see |Luj . 

Finally, a natural issue is the asymptotic estimates of the stable fibers as e ^ 0. As in 
Section m we use equation ()4.44p as the approximation of ()4.5p in the slow direction and 
keep the same notations as in (14. 43 p . Given any £ A"i> let = + (^s + ^u)('?c) £ 
where we recall that M.% is the center manifold of (I4.44p . ^* are independent 
of and x,,(^c)(i) be the solution on M.% such that x*(^c)(0) = C*. Let x^(t) satisfy 

i,(t) = e*^/6 + ( re(*-^)^/P.+ r e(*--)^/Pe.)F,^(i.(T),ec,T)dT. 

J+oo 

where _ 

(5, , t) = F^ {x + x,{C,),t + to)- Fl (x, (Cc) , t + to) . 
Therefore, (S* +x^,(.^c))(i) is the solution of the unperturbed fibre starting at the based 
point with height such that 

(5?, + x,)(0) =is + {I- Pu% + Ki^s + Pu)^*,e) ^ ct:J6, Cc, e) 

where /i* : Xf x x M ^ is the defining function of the center-stable manifold 
of (IQijl . 

Theorem 5.4. Fork = 2, assume (Al)-(A5), (B1)-(B5), (C1)-(C2), and A2 in ^Ml 
is nonempty. For = and the above given ^c, we have 
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where C depends on K,ai,a[,a2,r],r]' ,r,e^,, \£,c\xi, and \is\xi- Moreover, if k > 3 and 



XixYi 



<C'e{\gx\L{X,Y) + \^xGl\cO(XixR,L{X,Y)) + \9tGl\c0c^(^XixR,Y)) 

We will write instead of S^cy as we only consider the case = in this theorem. 
Proof. By (|5.4p and the definition of G 

/I 



(5.7) 

+ / e(*-^)(7+9.)(G(i,0,e„e)+5.i)dr. 

J+oo ^ ^ 

Since |D(j^jr)G|co = |D(2;,j,)G|co < r and ai < rj < rj + rj' < a2, 

(5.8) sup-e"(^+^')Vi(i)ln < „ , „ \y\n+v',^.,Y.- 

To estimate 12(^)5 we integrate by parts to obtain 

h{t) = e(J + egyY' ( - {g.,x{t) + G{x{t), 0, t, e)) + T e(*-^)(7+9.) 

((/a; + £'xG(x(Cc) + x,y{Q,T + to,e))x + dtG{x,0,Cc,T, e) 

+ (l)G(x(ec) + X, y{i,),T + to, e) - DG{x{ic),y{ic).T + to, e)) mc).y{ic)) 



dr 



where x, (x(^c)(''"), ?/(Cc)(''")), and their time derivatives in the above integral are all 
evaluated at r. Using the differential equations satisfied by x(^c) and x, it is easy to 
show 

(5.9) mc)\^',i,x + \K+v',i,x<G'. 

Since = 0, from ()4.5p and use (j4.50p to estimate its right side, we obtain for any r, 
G' 

(5 10) l2/(^'=)lv,i,v' <—\yiCc)\:!l>^i^Xi + ilakixx) + I^^G'^lco(XixR,y))|2;(Cc)|J,i,Xi 

<C''(lfi'x|L{Js:,y) + \DxGl\cO(X:t^xR,L{X,Y)) + \9tGl\c0ci{XixR,Y)) ^ 

where G' depends on K, a'^,r]' ,r, e*, ICcLyi. Using (j5.9p . the bound on DG, the estimates 
on ISIr;,!,^! and from Theorem 15.11 assumptions for /c = 2, it is straight 

forward to see Ihlrj+ri' ,ei, ,Yi — G'^ and along with (|5.8p . we have |y|r?+r?',e*,yi ^ G'e. 

In order to obtain a more careful estimate in terms of Qx and Gl when k > 3 and 
assuming the extra assumption on dtG, skipping t and e, we rewrite 

DxG{x{^c) + x,y{Cc)) = DxGl{x{Q + x) + [' DxyG{x{Q + x, syiCc))dsy{Cc). 

Jo 
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Similarly, rewrite G{x, 0, ^c) (as well as dtG and DG{x{^c) + x, . . .) — DG{x{S,c), ■ ■ •)) 

G{x,0,^c) = G{x{^c) + x,y{^c)) - G{x{Cc),y{^c) = C D^Gl{x + sx)dsx 

Jo 

+ 1^1^ D^G{xiCc) + six,S2y{Cc))dsids2[{x,0),{0,y{Cc)) 



Therefore, in the estimate of I2, each term either directly has a factor or Gl or 
indirectly from y{(,c) and (|4.5U|) which implies 

\y\ri+ri',e^,,Yi - G' ei\gx\L{X ,Y) + \^xGl\cO{XixR,L(X,Y)) + 1^* ^* Ic^Ci (Xi xM,y) ) ■ 

Using integral equations of x{t) and we have 

x{t) - x.{t) = ( ^* + £ Pc«)e(*^^)^/ (f{x, y, Cc, e) - ^^(5?*, Cc) + fyyjdr. 

We can write, skipping r + to and e, 

F{x,y,ic,e) - K{x*,ic) 
=F{x{ic) + X, y{i,) + y)- F{x{Q,y{i,)) - F{x, (^c) + x,,0)+ F{x, (C,, 0)) 

= / DF{x{(,c) + sx + {1 - s)x^,y{^,c) + sy)ds{x -x^,y) + 

Jo Jo Jo 

D'^F{six{^c) + (1 - si)x^{(,c) + S2X^, siy{(,c))dsids2(^{x{(,c) - x^{Cc),yi(.c)) , (^*,0)^ 

Combining the estimates on \x{(,c)-x*{S,c)\ri' ,i,Xi and \y{£,c)\ri' ,i,Yi (|4.5(J|) . | X* \ri,l,Xi 
from the standard theory (like in Theorem 15. ip . \y\rj+r]',€i,,Yi the above, we obtain 
the desired estimates on |5; — a;* and thus complete the proof. □ 

6. Normally elliptic singular perturbations to homoclinic solutions 

In this section, we will discuss the persistence of a homoclinic solution under nor- 
mally elliptic singular perturbations. We assume (Al)-(A5) for k = 2, (A5'), (B1)-(B5) 
in Section 4 and (C1)-(C2) after Theorem 14. 9i In this whole section, we assume 

(Dl) A generates a strongly continuous group on X and has finite dimension. 
(D2) There exist rj and rj' such that 

ai < 27] < r] < min{0, 02} , max{0, a'l} < rj' < 2r]' < a'2, 
ai < 1] + 1]' < a2 , ai+i]' <0, 

where ai,a2,a[,a2 are defined in (B5) and (C2). 
(D3) (j2.4|) has a homoclinic orbit Xh{t) such that |Ax;i(t)|xi is bounded and 

supe~°'^'^\xh{t)\xi < 00 , supe~''''2*\xh{t)\xi < 00. 

t>0 t<0 

(D4) There exists a invariant quantity H : Xi with DH € G\Xi,L{X, M)) 
such that 

H{0) = , DH{0) = 0. 
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(D5) At xo = Xft(O), 

DH{xo) ^ , dim{T.,„M^(~]T^,M^') = 1, 

where Mq and Aiff are the unperturbed unstable and center-stable manifolds 
of 0. 

Our goal is to study if ()2.ip has a homoclinic solution to when < e ^ 1 and 
how this problem are handled under normally elliptic singular perturbations. We 
will consider both the weakly dissipative and the conservative cases via a geometric 
approach based on invariant manifolds. For the former, a more analytic method based 
on the Lyapunov-Schmidt reduction may also work \Cli\ ISZ1| to give the persistent 
homoclinic solution, but the geometric method provide more information such as the 
transversality of the intersection of the stable and unstable manifolds. For the latter, 
we are not aware of such an analytic method even in similar regular perturbations, so 
we follow the geometric approach as in |SZ3j . 

Since / is independent of t when e = 0, in this section, we will write 

fo{x) = fix, 0, t, 0) goix) = g{x, 0, t, 0). 

We use Bp{p, S) to denote the ball in a space S of radius p centered at p which is often 
skipped if p = 0. We will also keep using Px and Py to denote the projections. 

In he following Subsection 16.11 a coordinate system around the unperturbed homo- 
clinic orbit. Subsection 16.21 is devoted to study the persistence of the homoclinic orbit 
under weakly dissipative perturbations and Subsection 16.31 is to study the conservative 
and autonomous case, i.e. /, g are assumed to be independent of t for all e > in 
Subsection 16.31 The example of the elastic pendulum will be revisited. 

6.1. Coordinates around the unperturbed homoclinic orbit. Locally near 0, 
we cut off the nonlinearity as in section 4 to obtain hcs , and thus the local integral 
manifolds. We also use A^^(to) to denote the global integral manifolds corresponding 
to the time extended by the flow from the local ones of systems ()2.ip and (12. 4p , where 
a = cs, u, cu, s,c , /3 = 0, e. The assumption H{0) = and the invariance of the fibers 
and H imply i^|xo = and thus T^Mq C ker{DH{x)) for any x € Mq- Moreover, 
by assumption (D3), TMq'* can be foliated into the disjoint union of invariant stable 
fibres which are with respect to the based point. Therefore, there exists a nonlinear 
projection € [Mq'^ , A4q) , which maps points on each fiber to their based point, 
such that 

riA,. = (0,0) , rug = /. 

The fiber invariance implies H = H o on Ai^f. So for any x G Mq and 6x G T^A^q*, 
using the assumption DH{0) = 0, we obtain 

DH{x)6x = DH{f\x))Df'{x)6x = DH{0)Df{x)6x = 0, 

which implies 

(6.1) T^M^q' C ker{DH{x)), Vx G M'q. 

Similar properties also hold for the unstable and center-unstable manifolds. 

In the enlarged phase space, we trivially extend the domain of H from Xi to Xi x Yi. 
Clearly, 

(6.2) H{0, 0) = , DH{0, 0) = 0. 
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To study the perturbation of the homochnic solution, we need to take a cross section. 
Let 

V = Axo + foixo). 

Since v € Xi C X, there exists a hyperplane S' C X that is transverse to v. Let 
E = (E' Pi Xi) X Yi, by using v G Xi, one can prove v and T, are transverse in Xi x Yi. 
Let QvtQ'v be the projections from Xi x Yi and X xY onto with kernel S and 
S' X y, respectively. We will identify the range of and Q'^, i.e. Mu, with M. Let 

-^0 = -^0 n(^o + ^) ' ^0 = -^o' 0(^0 + S) , X" = r,„A^^ , = T^.M^o"- 
Prom (j6.1|) . we have 

, Xi , C ker(L>ii'(xo)) Q S ^ H. 

We use Codimty(-Z^) to represent the codimension of a linear subspace Z in a Banach 
space W. On the one hand, since ker(Dff(0)) is a hyperplane, f € ker(Dff(0)), and 
f ^ S, we have Codims(n) = 1. On the other hand, (D4) implies X^ H^i = {0}- 
Moreover, dim-XJ* < oo implies Codim£(X;^ © Xi Yi) = 1. Therefore, 

n = Xf ©X^en. 

Let w e S be transversal to H such that DH{xo)uj = 1 and Q^), Qcs, Qu, Qy he projec- 
tions from S onto U},Xi ,Xi and Yi. Thus, 

S = span{uj} © n = span{a;} © Xf © X" © Yi. 

We will use coordinates 

= {Quj{p - xo),Qcs{p - xo),Qy{p - xq),Qu{p - xq)) 
={DH{xo){p - xo),Qcs{p - xq), Qy{p - xo),Qu{p - xq)) 
to represent any p G T, + xq. Locally, there exist 5 > and 

To : BsixT) ^ M X Xi , ^0 : i^^C^i) Xf, 

such that the graphs of Tq, ^'o are open subsets of ^Aff,^AQ, respectively. We extend 
Tq to Bs{X I ) X Yi trivially in y. 

To study the perturbed problem, let r be the cut-off radius defined in section 4, 
there exist ti > 0,t2 < such that 

r 



(6.4) |xi,2Ui < + ip^^i^i + \Dhu\co) + \Pu\il + \Dhcs\co)) ' 

where xi^2 = 2^/1(^1,2)- Recall that <l>(t,to)2:^ + Uj^) and <&^(t,2;) denote the flow maps 
with the terminal time t of (12. ip and (12. 4p . respectively. We first show that for any 
to S M, 7W^''(to) does intersect S near xq for e <C 1. 

Lemma 6.1. For any to £ I^i there exists a unique t' = t'{tQ,e) such that 
<^{to,to + t',x[,e) G xo + S, 

mto,to + t',x[,e) - xo\^^^y^ + \t' - til + l^tot'l < C'e, 

where 

x[ = x[{t', e) = PcsXi + hu{PcsXi,to + t', e) G n Xi 
and C" depends on constants in assumptions of this section. 
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Proof. The proof is obviously based on the Imphcit Function Theorem, however, we 
have to be rather careful due to the singular perturbation natural of the problem. We 
will use (9i$, to denote the differentiation with respect to terminal and initial time, 
respectively. For any to ^ let 

-f{t',e) = Q'^{^{to,to + t',x[,e)-xo) , -f{t' ,0) = Q'^{^\-t' , xi) - xq). 

Theorem 12.21 implies, for t' on any bounded interval, 

(6.5) |7(t',e)-7(i',0)| <C"e. 

To show the closeness of 7(-,e) and 7(-,0), using the definition of x'l, one can 
compute 

dt'j{t',e) - dtn(t',0) < Q'^D^{to,to + t',x'i,e)dtoK{x[,to + t',e) 



+ 



Q[, [D^to, to + t', x[,e)V,{to + t',x[)- D^^{-t', xi)Vo{xi: 



where V^it, x), Vb(x) represent the velocity field of (|2.ip and (j2.4p at (t, x), respectively. 
From Proposition 14.12] and Remark 14. 131 we have — xi|xi < Ce. Explicit computa- 
tions based on the forms of ([23]) and ([23]) imply \PxiV^ito + t' ,x[) -Vo{xi))\x < C'e. 
Applying Theorem 12.31 we obtain 



{D^ito, to + t', x[,e)V,{to + t',x[)- D$°(-t', xi)Voixi)) 



<C'e. 



From Theorem 14. 161 and Remark 14.13^ we have |i9to^ii(a;i, + t',e)\x < C'e and thus 



(6.6) 



<C'e. 



Q'^D^{to,to + t',x'i,e)dtohu{x'^,to + t',e) 

Therefore, we have proved 7(t', e) and 7(t', 0) are close for t' on bounded intervals. 
Since the system (12. 4p is autonomous when e = 0, 

7(ti,0) = 0, dtnituO) = -Q',di<l>\-ti,xi) = -Q> = -1. 

By implicit function theorem, there exists a unique t' = t'{to) such that 

<^{to,to + t',x'^,e) G xo + S \t'{to)-ti\ < C'e. 

Moreover, from Theorem 1 2 . 1 1 and the C^ smoothness of $(— t',xi) G Xi in t' which is 
due to the assumption Axh C Xi, it is easy to obtain 



\<P{to,to + t',x'^,e) 



<C'e. 



(6.7) 

Finally, note that 

dt,-f{t', e) =Q; (y,{to, $(to, to + t', x'l, e)) - D^{to, to + t', x'^,e)V,{to + t', x'^) 
+ D^{to,to + t',x'^,e)dt^,hu 



By Theorem 12.31 and (|6.6p . (|6.7p . we have \dtg^{t',e)\ < C'e, which implies the desired 
estimate on dt^t' . □ 

Next we consider the tangent space T(S P| A1^'*(to)) near ^(to,to+t'(to),x'i,e) based 
on Theorem 12.31 and more directly Proposition 14.15) and Lemma 13.31 Let E^{t,to,x) 
be the evolution operator defined in Theorem 12.31 for x G Xi, i.e. 

diE'{t,to,x) = {- + Dygoi<l>^it - to,x)))E'{t,to,x), E\to,to,x) = ly . 
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We notice that the operator E defined in Proposition 14.151 is only 0(e) away from E*" 
on any finite interval. Therefore, we have 

Lemma 6.2. Let t' = t'(to,e) be the one found in Lemma \ 5.1[ For any C > and 
small 6 > (independent of e) and {6.cs,(,y) G Bs{PcsXi, Xf'^) x Bctiyi), we have 

Moreover, if{6x,6y) G Xf* x Yi with \6x\xi + \Sy\Yi < 1, then 

($(to, to + t', ^cs + Cy + Ki^cs,^y,to + t',e),e))6x 
-D^^^{^^{-t',^,s + h'i{U))Sx ^ <C'e, 

and 

PxD/:^ {^{to,to + t', Ccs + Cy + hu{Ccs,Cy,to + t' , e) , e)) 6y £ Bc'e{X) n Bc>{Xi) 

[PyD^^ {Hto,tQ + t',^cs + ^y + K{Ls,iy,tQ + t' , e), e)) 
- ^(to, to + t'; U + hl{icu),e)) 6y 
where C depends on C and those constants in assumptions. 



<C'e. 



The smoothness of <I> and and the assumption dim— < oo, which implies 
the equivalence between | • and | • \x^, are used in the proof. Moreover, even though 
Proposition 14.151 is stated only for = 0, our assumption \S,y\Yi = 0(e) combined 
with the smoothness of in (Theorem I4.10p is sufficient to guarantee the above 
estimates. 

In the next lemma, we will write A^^* H {xq + S) locally near xq in the coordinate 
system {d,x^^ ,y,x'^). The main issues are the size of the coordinate chart of the 
manifold and the regularity estimates. 

Lemma 6.3. For any b > 0, there exist eo > 0, 6' > and 

T = (T^T") : Bb^ixT) x Bf,iYi) xRx (0,eo) ^ (M,Xi) 

such that xq + Graph{T{tQ,e)) is an open subset of Ai'^^{to) where 

Graph{T{to, e)) ^ {T"'(x-, y, to, e)^ + + ey + T"(x-, y, to, e)| 

x^^GS,K^r),yeB,(yi)}. 

Moreover, T are C"^ in x^^,y and satisfy 

(^•^) 1^ ~ ^0lci(Bi,'(^r)xB6{n),X5'xlR) + l^to^lco < C^'e, 

where C only depends on b and those constants in assumptions. 

For e = 0, we define T(x'^*, y, to, 0) = To(x'^^). Notice in the definition of Graph(T) 
we scale y to ey. This is to avoid the dependence on e of the domain where T is defined. 

Proof. The first step of the proof is to establish a correspondence between an open 
set of 7W^**(to), which is a hypersurface of 7W^*(to), near xq and a hypersurface of 
Mfito + t') near x'^. 
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Let w = Pcs{Axi + /o(xi)) and Xf C Xf such that Xf = Mu; Xf . For any 
61 > 0, define 

•^(a, Cs^^y: e) = ^{to,to + t' , P^sXi + aw + + 

+ hu{PcsXi + aw + i'^g,eiy,tQ + t',e),e) - xq. 

Here o G [—(5,5], ^ BsiXl'^) and G i3f,^(Yi), where 5 > sufficiently small but 
independent of e. From Lemma 16.21 we have 

^^■'^^ |-^('' ^) ~ ^)lci{[-<5,<5]xB4(Xf»)xB6^(yi),XixYi) - ^' 

where C depends on those constants in the assumptions. Lemma [6.1 [ implies Q^J-{Q, 0, 0, e) 
for all e G [0, eo). From the Implicit Function Theorem, we obtain that, when 5 and 
eo are sufficiently small there exists a : Bs{Xf^) x Bh-^{Yi) x [0, eq) — )> [—5, 6] such that 

(6.10) T{C,^y,e) ^ Ha{Cs.iy.e),i'cs,iy.^) G S, ^ {Cs.iy) € i^^l^f) x B,,{Y^). 

For e = 0, we identify a(Ccs')^y5 0) with 00(^^5)5 which satisfies 

Pc.xi+ao(ecs)^« + 4 + hliPcsXi + ao(eL)«^ + ^cs')) - G S n Xi. 

Moreover, by assumption (D2) and Theorem 14.21 o is in ^csj'^j/ > '^0 is in 
0(0, 0,e) = ao(0) = 0, and 

^^■^^) |a(-,-,e) - «o(-)|ci(Ba{jff'')xB6^(yi)x[o,eo),[-<5,<5]) ^C'e. 

Consequently, 

i^-^'^) l-^(T'e) - •^o(-)lcl(B5{Xf=)xBi,^(Yi)x[0,eo),Xixyi) - ^'^ 

where 

-^o(^L) =-^(ec.,ey,0) = .F(a(eL,ey,0),ecs',Cj;,0) = J^(ao(CL),^L,0,0). 

To obtain the estimate on dt^a, notice (j6.10p is equivalent to Q'^J- = 0. Differentiate 
it with respect to to and wrote S,cs for Pcs^i + aw + S^'^g, we note that 

\dtoQv^{a,Ccs^^y^(^)\ 

< g;(K(to, Ha, CCy, e) + ^o) - (1 + dtot')D^{to, to + t', Cos + eCy 

+ huiCcs, e^y, to + t', e),e)Ve{to + t' , Ccs + e^y + hu{Ccs, e^y, to + t',e)) + D^dt^^K) 



< 



C'e. 



g;(Vo(-F(a,cL,ey,o) + xo) - z?<&°(-ti,Ccs + /i°(ecs))Vo(ecs + 

Here we use (16. 9p . Lemma l6. 11 16. 2| Proposition 14. 16] Theorem 12.31 to obtain the above 
estimates. Finally, the term other than C'e in the above right side vanishes since the 
system is autonomous when e = 0. Therefore \dtQa\ < C'e which along with a similar 
procedure implies 

(6.13) \dt,T\x < C'e. 

We claim for any 6 > there exist 61, 6' > independent of e, such that the map 

{QcsJ", ^QyJ")-^ : Bb'iXl') X Bb{Yi) BsiX^n x Bb,{Yi) 

is well defined and its C^ norm is uniform in e. 
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To prove this, we need solve the equations 
(6.14) a.) QcsHi'cs,iy,^)=^''\ b.) -^QyT{^'^^,^y,e) = y. 

We first find a good approximation of this system of equations. By (|6.12p and Lemma 
16.21 one can compute 

QcsDc' T = QcsDt' Tq + eOi, QcsDcT = e02, 

where O = G L{X^^ xYi,Xi xYi) is bounded uniformly in e and E is the 

linear evolutionary operator defined in Lemma 16.21 at the base point 

e = Pcsxi + ao(CL)"' + + KiPcsXi + a^ii'Jw + CD- 
It implies that for fixed ^-iid e << 1, 

Since E and E^^ both have upper bounds independent of e and is C^, by (j6.15p 
and an Implicit Function Theorem argument, for any y G Bi^iYi), ^'^^ € Bs{Xf^) and 
e € [0,eo) b.) of (j6.14p has a unique solution CyiC'cs^V^^) ^^I'i l-^C — C". Along 
with (j6.15p . it implies that, as a mapping of 

Since QcsJ~o is independent of e and is locally invertible, one can use an Inverse Function 
Theorem argument again to prove there exist sufficiently small 6' > 0, eq > so that 
for {x'^^,y, e) G -Bfe'(0, Xi ) x Bi,{0, Yi) x [0, eo), there exists a unique icsi^'^" lU-, e) which 
is also in x'^'^ and y satisfying it along with iy{£,'csS, y, e) solve (j6.14p . Therefore, we 
proved the existence of {Qcs^, ^QyJ^)~^- The estimate on its norms follow directly 
from (j6.15p . 

For (x^^^/,^o,e) G Bb'iO,xT) x Bb{0,Yi) x M x [0,eo), let 

(6.16) + ey + T{x^',y, to, e) = /"((Qcs-F, iQJ,7•)"'(x^^ y), 

When e = 0, (|6.16p is becomes 

x- + To(x-) = 7-o((Q,,7-o)"'(x- 

Since J-" is C^, T is also C^. The estimates on DT follow in a straight forward manner 
by differentiating (|6.16p and using (j6.15p and (|6.13p . □ 

Finally, we present a similar coordinate representation of the stable manifold which 
is obtained in rather similar fashion. For any to € M, there exists a unique t" = t"{to, e) 
with \t" -t2\< C'e such that 

^{to,t" + to,x",e) Gxo + S, where x"(t",e) = PuX2 + hcs{PuX2,to + t" , e) 

with similar estimates as in Lemma l6.1i Moreover, D^(to, to + t" , e) satisfies similar 
estimates as in Lemma 16.21 except there is no D^^ terms. 

There exists 6 > sufficiently small but independent of e and 

= (^r'^,^-,^?/) : Bb(Tt) X M X [0,eo) ^ {R,xT,Yi) 
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such that Graph(^'(-, to, e)) is an open subset of 7W"(to) where 

GrapH^{-,to, e) ^ [x^ + to, e) + ^^(x", to, e) + M'^^(x", to, e)|x" G Bb(xl)}. 

Moreover, ^I*^ is in and satisfy 

(6.17) |^(-'*o,e) -^(^tcO)!^!^^^^^^^^) !^^^^.^^^^ < C'e , \dt,^\xxY < C'e, 
where C is independent of e. 

6.2. Persistence of homoclinic orbits under weakly dissipative perturbation. 

In this subsection, we assume additionally 
(A7) For i = 0, 1, 2, the following quantities have a uniform bound Co, 

{dt'D'f,dt'D'9) G C'iXi X n X R\Li{Xi x Yi,Xi x Yi)). 
In order to study the persistence of the homoclinic solution of (j2.ip . we first derive 
the Melnikov integral to measure the distance between 7M^*(to) and A^"(to). From 
the construction of the coordinate system, the intersection of Ai^^ito) and A^^(to) is 
equivalent to the following system: 

(6.18) = T"(x-, y, to, e) , x'^' = ^'^%x\ to, e) , ey = M'^(x", to, e), 

(6.19) d = T'^(x-,y,to,e) = ^'^(x",to,e). 

From ()6.8p . ()6.17p . and a contraction mapping argument, one can easily prove 

Lemma 6.4. There exists cq such that for every e € [0, eo) and to S M, there exist 
^cs _ x^^(to,e),x^ = x"(to,e),y = y{tQ,e), which are continuous in to and e, satisfying 
()6.18p . Moreover, 

(6.20) Ix'^^lxi + <C'e, 

where C depends on constants in assumptions and uniform in to and e. 
Let 

(6.21) P"(to, e) = (^''(x"(to, e), to, e), x"(to, e), x'='*(to, e), ey{to, e)) + xo, 

(6.22) ^^=^ (to, e) = (T'^(x^'^(to, e), y(to, e), to, e), x"(to, e), x^^(to, e),eyito, e)) + xo 

(x_(t), y_ (t)) ^ <I>(t, to, P"(to, e), e) (x+(t), y+(t)) ^ <I>(t, to, P'=^(to, e), e). 

From the coordinate system we constructed in the previous subsection, clearly, the 
center-stable and unstable manifolds intersect if 

(6.23) = P^' ^ ^'(x_,to,e) = T'^(x+,y+,to,e) ^ H{P^) = H{P'^), 
where (|6.3p is used. 

Melnikov method. From (j6.2p . we have 

H{P'-) = H{P^) - H{0) = r DHmt,to,P^e))dMt,to,P'',e)dt 

(6.24) = r Z5F(x_(t))(/(x_(t),y_(t),t,e)-/o(x_(t)))(it 

where the last equality follows from the fact that H is invariant under (|2.4p which 
implies, for any x G Xi, DH{x){Ax + fo{x)) = 0. 
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We claim for all t < to and max{a'^,0} < r]' < 
(6.25) \x^{t) - Xh{t - to)U^ + |y-(0|y, < C'ee"' 

In fact, from Lemma I6.2H6.4I and Theorem 12. 2^ this inequality is obvious for t G 
[t2 + to)*o]- For t < to + ^2, since (x_,y_) and Xh remain in a small neighborhood 
of the perturbed and unperturbed unstable manifold, respectively. From (I4.50p (and 
Remark I4.13P and the standard Lipschitz dependence of the unstable orbits on the 
initial base points in terms of the exponentially weighted norm (obtained from the 
uniform contraction mapping principle), (|6.25p follows. 

H{P^) = r DH{xh{t - to)) {Dyfy-{t) + edj)dt + 0{e^) 

where Dyf and d^f are both evaluated at {xh{t — to),0, t,0) and the smoothness 
of H and DH{0) = are used to guarantee the convergence of the above integral. To 
estimate the y-(t) term, we write the variation of parameter formula, 

y_(t) = e(*-*«)^y(to)+ / e^'~^'>i g{x.,y.,T,e)dT. 

J to 

On the one hand, integrating by parts and using (|6.25p . we obtain 

r DH{xh{t - to))Dyf{xh{t - to),0,t,0)e^'-'^^iy{to,e)dt = 0{e^). 

On the other hand, we compute by changing the integration order, integrating by parts, 
and using the exponential bounds of the orbits in t, 

(6.26y" " DH{xh{t - to))Dyf{xh{t - to), 0, 0) ( ^ e^'~^'^ig{x. (r), y_(r), r, e)dr) dt 

(6;S7-> r (f DH{xh)Dyf{xh,0,t,0)e^'~^'>idt)g{x.,y.,T,e)dT 

{6^8^e [" DH{xh)Dyf{xh,0,T,0)J~^g{x^,y^,T,e)dT 

( y -{DHixh)Dyf{xh,0,t,0))j-^e^'--Hd?jg{x_,y_,T,e)dT 
It's easy to see from (j6.25p 

(6.29) /" DH{xh)Dyf{xh,0,t,0)J-^g{x_,y_,T,e)dT 
J —00 

= - r DH{xh{t - to))Dyf{xh{t - to), 0, t, 0)J-^go{xh{t - to)))dt + 0{e^). 
J —00 

Again, integrating by parts on e^*~^^"^ and using the assumption Axh E ^1, one may 
compute 

( y -{DHixh)Dyf{xh,0,t,0))J~^e^'~-^idtjg{x_,y_,T,e)dT = O{e). 
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Summarizing all the estimates, we obtain 

(6.30) H{P'') = ej uj{t,to)dt + 0{e'^), 

J —oo 

where 

(6.31) u{t, to) = DH{xhit)){dJ - DyfJ-^g) {xh{t), 0, t + to, 0)) . 
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Even though (^{t,to,P'^'^{tQ,e),e) does not necessarily stay in a small neighborhood 
of the origin for all t » 0, we will still obtain a similar approximate for H{P'^^). In 
fact, by the same argument leading to (|6.25|) . one can show that an inequality similar 
to (|6.25p holds for t € [to, 2^] as long as (x+(t), y+(t)) stay in the r neighborhood of 
for all t E [to + ii , r] , where r is the cut-off radius in the construction of the center- 
stable manifold. Therefore, let a = < and Ti = aloge > 0, we have, for any 
t G [to,Ti+to], 

(6.32) 



{x+{t),y+{t)) - Xh{t - to) 



where C is independent of t and e. In particular, 

(6.33) \ix+,y+){T, + to)-xoiTi)\x,>cY, < \xhiTi)\x, < C'e"'^^ 
and 

(6.34) aai = l + ar]' >^ , 1+ 2ar]' > 0. 
Since H{0) = 0, we can compute 

(6.35) H{P''') - e / uj{t,to)dt= H{x+{Ti + to)) - e / u}{t,to)dt 

J +00 J +00 

.0 

+H{P^')-H{x+{Ti+to))-e / u{t,to)dt. 
Using (j6.33p and a similar procedure as in the approximation of H{P'^), we obtain 



(6.36) 
(6.37) 



H{P'') - H{x+{Ti + to)) - e r uj{t,to)dt < C'e^^^"'^', 

JTi+to 



uj(t, to)dt 



+ 00 



H{x+{Ti+to)) < \D^H\co\x+{Ti+to)\' < C'e^+^-^'. 

Therefore 

(6.38) F(P") - H{P^%e)) = eM(to) + 0(e2+2'^'''), M(to) = 



+00 



uj{t, to)dt. 



By Lemma 16.41 (j6.23p . and an Intermediate Value Theorem argument, we obtain 

Lemma 6.5. Suppose M(to) has a simple zero at some to, then there exists eo such 
that for each e € [0,eo), there exists t* such that 7W^*(t*) and 7W^(t*) intersects near 
xo = Xh{0). 
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Remark 6.6. Compare the Melnikov functional obtained in the above with the one 
under regular perturbations, we observe an extra term —DyfJ~^gQ. In fact, it is easy 
to see where this term comes from the coordinate change as in ()2.5p . Let 

XI = x, yi = y + eJ~^go{x). 

Then the y equation takes the form of (j2.6p where gi = 0(e). From Proposition 
\4-12 and \4-15 and their remarks, it is easy to prove that the contribution from the y 



equation to the invariant manifolds are of order O(e^) and it does not appear in the 
leading Melnikov functional. The x equation now takes the form of 

xi = Axi + /(xi, yi - eJ^^go{xi), t, e). 

Ignoring yi , the Melnikov functional of this regularly perturbed equation is exactly the 
one obtained in Lemma \6. 5\ The only reason we did not take this approach is that thus 
transformation reduces the smoothness of the system by 1 order which would require 
k > 3 in the assumptions. 

Homoclinic solution. Lemma 16.51 gives a condition for nonempty intersection of 
center-stable and unstable manifold. This intersection means the existence of a solution 
which converges to the steady solution as t — >■ — oo. As t increases and t < a log e + to, 
based on the stable foliation in the center-stable manifold, this solution will approach a 
neighborhood of the steady state inside the center manifold. In order to find conditions 
for this solution to converges to the steady state as t — )> oo, in this subsection, we 
focus on the case when the unperturbed center manifold is at least neutral and the 
perturbation is weakly dissipative so that the perturbed center manifold is weakly 
stable. In this case, the size of the basin of attraction of on the center manifold 
is the key issue. We illustrate how the method in the regular perturbation cases can 
be adapted here under certain assumptions, which are not optimal as we are only 
giving an illustration. In order to specify the assumptions, we first look at the Taylor's 
expansions of / and g, 

f{x,y,t,e) = f{x,y,t,0) +ef2{x,y,t,e) = fxX + fyy + fi{x,y) + ef2{x,y,t,e), 

gix, y, t, e) = g{x, y, t, 0) + eg2{x, y, t, e) = g^x + gyy + gi{x, y) + eg2{x, y, t, e), 

where fx^y = Dx^y f {0,0, t,0) and gx,y = Dx^yg{0,0,t,0) which are independent of t. 
Let Pc^su be linear projections from Xi onto x^'^" which are invariant under e*^'^+'^''\ 
where Xi = X^ © Xf". For any x S Xi, we denote Xc = PcXi and Xsu = Psu^. In 
addition, we assume 

(El) dimXf < +00 and {f,g) are in {x,y) with upper bound uniform in t. 
(E2) Let A{e) = A + fx - e and = ^ + gy - e. We further assume a[ < in 
(C2), which implies 

\e'^^'^\LiXf) + \e'^y\LiY^)<Ke-'' for t > 0, G 
(E3) For {x„x,,u,y,t,e) G X^ x Xf« xY^xRx [0,eo), 
Pcfy = 0, Dl^y^PMO,0,0)=0, 

Pcf2{xc,Xsu,y,t,e) = -Xc + eBo{e){xc,Xsu,y) + Bi{xc,Xsu,y,t,e), 
Bq{€) is a bounded linear operator acting on (xc, Xsu,y), 
Bi{0, 0, 0,t,e) = , DBi{0, 0, 0, t, e) = 0. 
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(E4) For {xc, xsu, y, t, e) G x Xf« x x M x [0, eo), 
5x = , gi (0, 0, 0) = , D^^^y^gi (0, 0, 0) = 0, 
g2{xc, Xsu,y, t, e) = -y + eB2{e){xc, Xsu, y) + B-sixc, Xsu, y, t, e), 
i?2(e) is a bounded linear operator acting on {xc, Xsu,y), 
B-siO, 0, 0, t, e) = , DB^iO, 0, 0, t, e) = 0. 

Remark 6.7. It looks that the above assumptions are too restrictive. However, one 
should first try to 'diagonalize' the linear part to remove fy and gx- (As a separate 
topic, we will discuss this transformation in the Appendix.) With this 'diagonalized' 
linear part, one is in a position to carry out a normal form transformation to elim- 
inate some quadratic terms. Assumption (E) should he considered for the form after 
performing a normal form transformation. 

For sufficiently small r, from Theorem 14.11^ for e G [0, eo) and to ^ there exists a 
local center manifold 7V4^(to) which contains as an open subset the graph of 

hsu = {^u.'^s) : Br{Xl) X Br{Y{) X [0,eo) ^ x 

where hsu{-,-,^) is understood as independent of y. Moreover, from Section [U hgu is 
uniformly bounded in in {xc,y) and Proposition 14.121 and Remark 14. 141 implv 

/is«(0, 0,e) = , \Dhsu{xc,y,e)\ < C'(e+ \xc\ + \y\) 

(6.39) \hsuixc,y,e)\ < C"(e+ + \y\)i\xc + \y\) < r, 

{x, + y + hsu{xc,y,e)} = M'^.ito) n {BriXD x Br{X^) x B^iXf) x i?,(Yi)), 

where C depends on those constants in assumptions. Here the assumptions that X" 
and X^ are finite dimensional are used. On the center manifold, the flow is reduced to 
the Xc and y direction only, where the solutions are given in the form of 

Xc{t) =e(*-**)^(^)x,(t.) +^*e(*-^)^(^)7(x„2/,T,e)dT, 



y{t) =e(*-**)4^y(t.) + [\'^'-^^^g{x,,y,T,e)dT, 



(6.40) 



where _ 

f{xc,y,t,e) = Pcifi + e/2)(xc + hsu{xc,y,e),y,t,e), 

g{xc, y, t, e) = {gi + eg2){xc + hsu{xc, y, e),y, t, e). 

From assumptions (El) — (E4) and (|6.39|) . {f,g) satisfies 

7(0, 0, t, e) = , g(0, 0, t, e) = , \Df\ + \Dg\ < C'{e^ + \xc\^ + \y\^). 

Suppose we have solution {xc{t),y{t)) such that |xc(t^)| + |y(t*)| < Se^ and |xc(t)| + 
\y{'t)\ ^ (2 + K)6e2 for t G [t*,^'] for some T' . By using Gronwall's inequality, we 
obtain from ()6.40p 

e'^'-'*\\x,{t)\ + \y{t)\) < K5eh^'^'<'-'*1 

Since C is independent of t and 6, by taking 6 < , we can extend T' to +oo and 
derive 

\x,{t)\ + \y{t)\<K6eh--2<'-'*\ 
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Therefore the basin of attraction on the center manifold, and thus on the center-stable 
manifold as well, contains the graph of hgu over the ball B i {Xf x Yi). 

Recall that (x+,y+) denote the solution with initial time t* we obtained in Lemma 
16.51 which satisfies 

In particular, for sufficiently small 6 and by choosing t-^ = Ti + t* , (j6.33p implies for 

t>Ti + t* 



\PcX+it)\ + < K5e2e 



■\e{t-Ti-t') 



— > as t ^ +00. 



Theorem 6.8. Assume (A1)—-(A4), (A5'), (A7)fork = 2, (B1)—(B5), (C1)—(C2), 
(Dl) — (D5) and (El) — (E4-)- Suppose the Melnikov function M(to) has simple zero 
points, then there exists eo > such that for any e G [0, eo), (|2.ip has homoclinic 
solutions to the origin. 

Elastic Pendulum revisited. Finally, we would like to revisit ()1.4p . We let y = eu 

and y = ui — e'^^u. We rewrite ()1.4p as a first order system. 



(6.41) 



We assume 
(PI) 



(1 + en)2 

xi = —g{l + eu) sinx — 2e'yxi + eFi(x, eu, t, e) 



u = -ui — eju 
e 



Ul 



1 X 

-u - ejui + / g 
e (1 + euy 



+ £ 7 u + g cos X + eF2 {x, eu, t,e). 



7>0 , Fi (tt, •,*,£) = , dtF2iTT,-,t,e) = 



where the assumptions on F are for simplicity. By implicit function theorem, there 
exists a locally unique steady state (vr, 0, n*^, lif ) which satisfies (u^,nf) = (0(e),0(e^)) 
and 



(6.42) ul - e^n" = , -u^ - e'^-ful + e^7^n' - eg + e^F2(7r, eu',t, e) = 0. 

Let {x,v,vi) = {x — t: ,u — u'^ ,u — u\) to translate the steady state to 0, ()6.4ip becomes 



(6.43) 



X =Xl + 



1 



l)xi 



1 + en*^ + evY 

xi =g sin x + eg{v + u^) sin x — 2e^xi + eFi (x + vr, eu^ + ev,t,e) 



V =-vi — e-fv 
e 

1 x^ 

v^=--v- e-fvi + — -3 

e (1 + eu^ + ev)-^ 



+ e^7^f + g{l — cos x) 



+ e(F2(x + vr, eii^ + ev,t,e) - F2{tt, eu^ ,t, e)). 
We rewrite the right hand side of last equation in (j6.43p as 

--V - ejvi + eDa;F2{7r,eu^,t,e)x + g{x,xi,v,vi,t,e) 
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where in view of (|6.42|) 

^2 ^ ^ 

g{x,xi,v,vi,t,e) =-— j- + g{l -cosx) + e^j'^v + e(F2{x + 7r,eu'' 

(1 + eu^ + evy ^ 

+ ev,t,e) - F2(7r,eti^,t,e) - DxF2{TT,eu'' ,t,e)x) . 

Clearly, (16.430 does not satisfy assumptions (E1)-(E4) due to the presence of the 
linear term €DxF2{'Tr,eu'^,t,e)x. We will eliminate this term by a linear coordinate 
transformation. A more general procedure of this type of transformation can be found 
in Appendix [3 To simply our notations, let 

\g Oj ' V-i -ejj e ^ ' VL>^F2(7r,eu%i,e) Oj 

By implicit function theorem there exist Li £ L{M?,'M?) with \Li\ < C'e^, where C 
depends on constants in assumptions, such that 

JLi - e^jLi - eLiMi + e^Mg = 0, 

which implies 

// OVVMi 0\(I 0\ /Ml 0\ 
\Li I J [eM3 M2j\Li I J \0 Ah J ' 

When e = 0, is a hyperbolic fixed point of the first two equations of (j6.43p . thus, all 
assumptions in (E3) for / are automatically satisfied. It remains only straight forward 
verifications that (E1)-(E4) are satisfied. Therefore, one can apply Theorem 16.81 

6.3. Persistence of homoclinic orbit under conservative perturbation. In this 
subsection, in addition to those assumptions given at the beginning of the section, we 
further assume 

(D7) Xf is finite dimensional. Moreover, 

dtf{x, y,t,e)=0 , dtg{x, y,t,e)=0, dim(r,„Xg f| T^o^Ig") = 1. 

(D8) There exists a family of invariant quantities {H{-,e) for (12. ID which, in terms 
of the Taylor expansion in u = |, takes the form 

H{x,eu,e) = HQ{x,e) + Hi{x,e)u + H2{x,e){u,u) + H^{x,u,e), , Hq{x,{)) = H{x) 

H, G C^-\Xi X M, Li{Yi,M.)) for i = 0, 1, 2, € C^{Xi x Fi x M, M). 

Here with a slight abuse of the notation, we still denote the unperturbed in- 
variant functional by H{x). Moreover, we assume there exist cq, C2 > 0, ci > 0, 
such that for any G X^ and {x,u) G Br{Xi) x BiiXi) 

Ho{0,e) = 0, DHo{0,e) = 0, D^HoiO,0){Ce,^c) > co\Ce\\ 
Hi{0,e) = 0, \DHi{0,e)\<ci , \H2iO,e)iu,u)\ > C2\uf, 

H3{x,0,e)=0, DuHs{x,0,e) = , \D^Hs{x,u,e)\ < Coe , a ^ coCa - c? > 0. 

Under these assumptions, one may still compute the Melnikov functional, but mostly 
it turns out to be identically zero. Our goal is find the intersection of the center-stable 
and the center-unstable manifolds. 
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We first refine the coordinates on the cross section S defined in Subsection 16. 1[ Let 

Much as in Subsection 16. H we have 

n ^ S n ( keT{DH{xo) © Yi) = © © © X^. 

Recall we took u € S\n with DH{xq)uj = 1. Let Quj,y,s,u,c he the projections on S 
given by S = Ma;©n and the above decomposition. For any p € S + xq, its coordinates 
can be written as 

(d, ={Quj,Qs,Qc,Qy,Qu){P - Xq) = {DH{xo),Qs,Qc,Qy,Qu)iP - Xq). 

Similar to that the center-stable and unstable manifolds in S can be written as graphs 
of T and ^ as given in Lemma 16.31 for any 6 > 0, there exist r > and Ti(-,e) : 
Br(xT) X M^i) ^XfxR and ^i(-,e) : Br(xl) Xf" x x M such that 

{ (Tf + TJ)(x^ x", y, e) + + x" + ey} C Mf, 

and Ti,^i satisfy similar properties as T,^ in (|6.8p . and (|6.17p . 

To find the intersection of A^^* and Aif^, we first try to match all coordinates except 
the d direction. Given (x*^,?/) € Br{Xi) x Bi,{Yi), by using the Contraction Mapping 
Theorem and Lemma [631 we obtain a unique pair x^'^{x'^,y,e) such that 

(6.44) T"(x^ x^(x^ y, e), y, e) = x"(x^ y, e) x^(x^ y, e) = TUx', x"(x^ y, e),y, e) 
and they satisfy, for some C' independent of e. 

x'*'"(0,-,0) =0 , i?^cx''"(0,0,0) =0 , |x"'"(-,-,e) -x"'"(-,-,0)bi < C'e. 

Among the above points on M'^^, next we identify the one on M^. More precisely, 
substituting x*(x'^,y,e) into T,^i and using Contraction Mapping Theorem, Lemma 
16.31 and an inequality for similar to (16.17p . we obtain a unique pair (x'^(e), y(e)) 
such that 

(x'=(e),y(e)) = {^l,-^^l){x'ix'{e),y{e),e),e), \x'^ie)\xi + |y(e)|y, < C'e 
which implies 

T(x'=(e), x^(x^(e), y{e),e),y{e), e) + x'{e) + x'{x'{e),y{e), e) + ey{e) 
(x^(x^(e), y(6), 6), 6) + x'{x'{e),y{e),e) G Ml 
Similarly, there exist (x^(e), ?/i(e)) = 0(e) satisfying 

Ti(x^(6),x"(x5(6),yi(6),6),yi(6),6)+x5(e) + x"(x5(6),yi(e),e) + eyi(e) 
=^ (x^ (x^ (e) , 2/1 (e) , e) , e) + x" (x^ (e) , yi (6) , e) € 
Let, for T e [0, 1] 

qir) = {qc{r),qy{r)) ^ {I - r){x\e),y{e)) + r{x\{e),y,{e)) 

p'^ir) = Ti{q{T),x^{q{T),e),e)+q,{T)+x^{q{T),e)+eqy{T) ^ x^{r) + eqy^r) G MT 
ptir) = T{q{T),x%q{T),,),,)+q,{T)+x'{q{T),e)+eqy{T) ^ xf (r) + eg,(r) G A^f . 
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We will show there exists tq € [0,1] such that p"(ro) = pf(ro), which is equiva- 
lent to H{x^{tq)) = H{xI{tq)) as we have matched all other coordinates in (j6.44|) . 
Since DH{xo)\j^ ^ and |y(e)y^, |yi(e)|yi = 0(e), by assumption (D8), it is clear 
that H{x^{to)) = H{x^{to)), and thus p"(to) = P^tq), if and only if if(p"(ro),e) = 
HiptiTo),e). Let 

h{T) = H{pt{T),e)-H{p^{T),e). 
To analyze /i(r), note (D8) implies ^^^(0,0) = , DH,{0,0) = 0. For any 

z = ^^ + eu + hsuiCc, eu, e) = x + eu e {Br{Xi) x B,b{Yi)) n M",, 
by using ()6.39p . we have 

\hsu\xi < C'{\Cc\xi{e + \^,c\xi) + e^|'u|yi(l + \u\yi)) 

Hoix,e) > (| -C"(6+ |eek))|Cc|i, -C'e^lnll.^, 

where C depends on the constants in the assumptions and b. Moreover, 

|^i(x,e)n| < ci{l + C'ie+\Cc\x^))\Cc\xMYi+C'e^\u\l^, 

H-2{x,e){u,u) > {c2 - C'{e^\u\Yj + \e.c\xi))\u\Y-,, \H3{x,u,e)\ <Coe\u\Y^. 

By the last inequality in (D8) and choosing sufficiently small r, there exists c* > 
such that 

H{z,e)>c,{\Cc\^ + \u\l^J. 

It implies that H{-,€) > in M'^f] {Br{Xi) x Bi,^{Yi)) except at 0, with quadratic 
lower bound (after the scaling y = eu). Actually, it also implies the origin is stable 
both in forward and backward time on the center manifold. Consequently, A4f are 
unique, where a = c, cu, cs. From the invariance of H{-, e), H{-, e) > in . 

From the definition of h, it is clear h{l) > > h{0). By the Intermediate Value 
Theorem, /i(ro) = for some tq G [0, 1]. 

Theorem 6.9. Assume (Al)—(A5), (A5')fork = 2, (B1)—(B5), (C1)—(C2) and 
(Dl ) — (D8). There exists eo > such that for any e G [0, eo), the center-stable manifold 
and center-unstable manifold of (|2.ip has nonempty intersection. 



The intersection of the center-stable and center-unstable manifold is generically 
transversal and forms a high dimensional tube homoclinic to the center manifold. See 
|SZ3j for more discussion in the regular perturbation case. 

Elastic pendulum revisited. Assume the elastic pendulum system (|1.4p is conser- 
vative, i.e. 7 = and the perturbation {eFi,eF2) comes from a small perturbation 
eG{x, y, e) to the potential energy. System (jl.4p becomes 

Xl 



(6.45) 



+ 



2 



Xl = -g{l + y) sin X - eD^G{x, y, e) 



ey = yi eyi = -y+ , + e^g cos x - e^DyG{x, y, e) 



{i + yf 

Its energy is given by the sum of the kinetic energy, gravitational, elastic, and pertur- 
bational energy 

^ = ^(1 + yf4 + ^ + ^ - 5(1 + y) cosx + eG{x, y, e). 
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Prom the Implicit Function Theorem, for each e, there exist a unique fixed point 
{x%0,u^,0) with {x%y^) = (vr + 0{e),ge^ + 0{e^)) such that 

(6.46) g{l + y')smx' + eD^G{x',y',e) = - e^g cos x + e^DyG{x\y^ ,e) =0, 

Let X = X — x"" ,y = y — y'' , we can rewrite ()6.45p as 

^ _ Xl 
''-{l + ye + y)2 

Xl = - g{l + y" + y) sin (x*^ + x) - eDxG{x'^ + x, y'' + y, e) 
^ 1 

y =-yi 

e 

1 1 ex^ 

yi = - -y- -y" + tT, — + ^gcos (x^ + x) - e^DyG{x^ + x,y^ + y,e), 

whose invariant energy takes the form 



— G(x'^, y^, e) j — ^[^(l + + ev) cos (x*^ + x) — (1 + y*^) cos x* 
where f = | , ui = ^ = y. Its Taylor's expansion yields 

^0 =^T7T^Sw - ai.'^ + ( + ^) - + ^('^(^' + y% e) - G(x^ e)) 

2(1 + y"^)^ 

^1 =( - n ^""^Ns + - - e^cos (x^ + x) + e^DyG{x' + x,y',e),0), 
H, = (l + YDlGix^ + x,y"'^^^^^ Oj. 

One can use (16.460 to verify the above Hi, where i = 0, 1, 2, satisfy assumption (D8). 
Therefore, for e <C 1, the center stable manifold and center-unstable manifold of (16.450 
intersect near the unperturbed homoclinic orbit Xh{t), which generically form a 2- 
parameter family of solutions homoclinic to a small neighborhood of the fixed point 
on the center manifold which is foliated by periodic orbits corresponding to small 
amplitude fast oscillations. 



7. Appendix 

In the appendix, we outline a procedure to block diagonalize the linearization of 
(|2.ip as a steady state via a linear transformation. And we will discuss two cases, 
namely, ^4 is a bounded linear operator on X and A is generator of a semigroup on X. 
This lays the foundation for further normal form transformations to eliminate some 
nonlinear terms. Assume 

(B) For (t,e) € M x [0,eo), 



(/,5)(0,0,t,e) = , dt{Df,Dg){0,0,t,e) =0. 
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The linearization of (|2.1|) at is given by the operator 



A, 



'A + D^f{Q,e) Dyf{0,e) 
D,g{0,e) j + DygiO,e)J- 

We look for linear operators L| € L{Xi, Y) and L| E L{Yi, X) such that their graphs 
are invariant under Ae- For simplicity, we write Dx,y{f, g){0, e) as Dx,y{f, g). Therefore, 
(L|,L|) should satisfy the following system: 

{J + €Dyg)L\-€L\{A + Dxf + DyfL\) + €Dxg = {), 

LlieDxgLl + J + eDyg) - e{A + Dxf)Ll - eDyf = 0. 

When A is bounded, treating the above as the equation for zero points of a mapping 
from L{X, Yi) x L{Y, X) x M to L{X, Y) x L{Yi,X), simply from the Implicit Function 
Theorem, we obtain 

Lemma 7.1. Assume A € L{X,X), (A2) and (A3) for k = 1 and (B). There exists 
eo > such that, for any e G [0, eo); (|7.ip has a unique pair solution of 0{e): 

{Ll,Ll)£L{X,Yi)xL{Y,X). 
Let = as the new variables. System (I2.ip becomes 

( S: = {A + Dxf + DyfL{)x + F{i, y, t, e) 

(7-2) S . J 

\^y = {- + Dyg + D^gLDy + y, t, e). 

Using the following embedding of spaces 

L{X,Y^) C L{X,Y) , L{Y,Yi) C L{Y,Y) , L{Y,Yi) C L{Yi,Yi), 

it is easy to show F and G satisfy propertie (A3) and thus all results in previous 
sections still hold for the new system (j7.2p . 

In the case when A is unbounded, we use an integral equation appraoch to solve 
(j7.ip under the assumption 

(F) There exist closed subspaces X*''^''^ of X such that X = X-' ® X" ® X"" and 
A + Dxf{0,e) is invariant on X''"''^. Let A^'"'' = + L>^/(0, e)) We 
further assume A^ is a bounded linear operator on X'^ and there exist uog < 
and uju > such that 

|e*^'| < Ke^^* for t>0 , |e*^"| < i^e"^"* for t < 0. 

To find Lf, we consider 

/■O 

(7.3) G"(L)x"= / e"{eLDyfL'' -eDxg-eDygL'')e-'^^\''dt, x"" eX"" 

J +00 

(7.4) G"(L)x" = e( J + eDyg)-\LDyfL^ - D^g + L'A^)x^ x'^ e X'^ 

(7.5) G^(L)x" = / e"{eLDyfL' - eD^g - eDygL')e-'^^' x'dt, x' G X^ 

J —oo 

where L'*'^'* = LP^^cs aiid -Pu,c,s are the projections defined by the decomposition. Let 
G{L) = G''{L)Pu + G%L)Pc + G'{L)Ps. 
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Lemma 7.2. // (j7.7|) below is satisfied, then G is a contraction from a bounded ball 
in L(X, Y) to itself and its unique fixed point L satisfies the first equation of ()7.ip and 

mL(X,Y) - ^' ' mL(Xi,Yi) - 
Proof. In this proof, we use the equivalent norm on X defined by 

||x|| =max{|P„x|, Pcx\, \Psx\}, 
which induces equivalent norms || • || of operators. In particular, 

11-^11 = \^^'\l{X^,Y) + \L^\l{X'=,Y) + \L^\l{x^,y)- 

Let 

B = {l|||L"|| < l}, Co = max{||A^|l, \\Dyf{0,e)\\, \\D^g{0,e)\\, \\Dyg{0,e)\\} . 
G satisfies the following estimates on B{p) 

(7.6) ||G(L)|| < 5GoK(-^ + r^ + 

\\G{Li) - G{L2)\\ < mCo(j^ + A + ll^i " ^2||. 

^I^^sl \^u\ ' 

If uju.s satisfy 

(7.7) J_ + J_ + 2e|J-i|<-^, 

\^s\ V>u\ 5a Co 

then clearly G defines a contraction mapping on B. Therefore, G has a fixed point L. 
Integrating ()7.3p by parts 

=r^{eLDyfV - eD^g - eDygL")x'' 

+ / J~^e"{eLDyfL''-eD^g-eDygL^)e-'^^\A^x''dt 
shows G L{Xf,Yi) and 

(7.8) JL" = eLDyfL"" - eD^g - eDygL^ + eL"^" 
which immediately implies 

mL{X]',Yi) - 

Identity (j7.8p and the similar one for along with the definition of G'^ implies L 
satisfies □ 

To solve for Lli consider L € L{Yi,Xi) with L"'^'* = Pu^c,sL and 

e / e-^*^" {Dyf - L^D.gL - L^Dyg) e'^ dt 

J +00 

'G(L) = <! e{Dyf - L'D^gL + A'L) (J + eDyg)-^ 

e f e-^*^^ (Dj,/ - L'D^gL - L'Dyg)e'^dt. 

J — oo 

Lemma 7.3. If 0Ju,s satisfy ()7.7p . i/iere exists a unique L € L(yi,Xi) sitc/i t/iat 
G{L) = L. Moreover, L € L{Y,X) with estimates \L\i(^Yi,Xi) ^ C' , |-Z>|L(y,x) ^ C^'e. 
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Proof. The proof of the lemma follows from the same procedure on the set 

Bi = |l \L^\l{Yi,X]^) + \L'^\l{Yi,X'^) + \L'^\L(Yi,Xf) < l|- 

□ 
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